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Department of Theoretical Physics, Bolestaw Bierut University, Wroclaw 


(Received December 1, 1956) 


In this paper the behaviour of the higher factorial moments of electron-photon 
cascades at large depths is investigated in approximation “А”. The differential cross 
sections satisfying the conditions of homogeneity are positive, and, in addition, are 
arbitrary functions. It is assumed that the Bremsstrahlung cross section is greater than 
the pair formation cross section. 

Тһе main results of this paper are: ` 
1) The main factor determining the development of the cascades at large depths is the 
product of the total number of particles and the threshold energy. Systems with a smaller 
value of this variable are privileged. d 
2) Photons, as primary particles, are *stronger'- than electrons — generalization of the 
conclusions of Stachowiak (1956) for moments of any order. 

3) Privileged configurations of electrons and photons were obtained. It turns out that 
configurations containing an even number of electrons different from zero are privileged. 
The weakest configuration for М particles is N photons and zero electrons when М > 1. 


$ 1. Introduction 


The object of this paper is to investigate the asymptotic relations satisfied by 

the higher factorial nioments of electron-photon cascades for large depths. The problem 
_ is studied in the Approximation “А“, і. e. the ionization losses being neglected. We 
assume that the cross sections satisfy the following conditions: 1) The differential 
cross section for the emission of a Bremsstrahlung photon per unit length has the form 
ил (ЕЈЕ,) dE] E, 
where E, is the electron energy before the emission of a quantum 

E is the electron energy after the emission of a quantum 


2) Тһе differential cross section for electron-positron pair formation in a unit length 
has the form w*(E/ Ej) d EJ E, 
where E, is the photon energy 

Ё (305) 
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E is the energy of one of the electrons (i. e. an electron or positron) arising 
as a result of the reaction. 


3) ОЎ (el и? (е) > 0 for O<e< 1 

1 1 
4) ] wt (=) d & = 0 — Ја (0 4: = ә 
0 0 

We adopt the following notation: 

А) 4,» (21... End Expt Entms X) 481... de, ,,, is the probability of finding at de- 
33 €; + deg) and m pho- 
tons in the infinitesimal interval of energy (e,; ¢,+ de, ), і = 1, 2,... n; Е = n + 1,... 
n-+ m and any number of electrons and photons of any energy, if the primary particle 
is, for i = 1, an electron of unit energy and for i = 2, a photon 2 unit energy. The 


pth x n ТАЕ, in the infinitesimal Е of energy (ғ; 


quantity @ у= bears the name "product densities". 

B) P'(n, m|e, x) is the probability of finding at depth x n electrons and m photons 
of energy greater than e and any number of electrons and photons of energy smaller 
than e, if the primary particle is, for 1 = 1, an electron of unit energy, and, for i = 2, 
a photon of unit energy. 


C) со со 
КИТТЕК = > У? H E P'(k, те, x) 


k=l п=т 


is the factorial moment. 


D ; 
) E итше х) 2, 2] P! (n, m | в, x) ии" 


п=0 т=0 


is the generating function of the probabilities P'(n, т | ғ, х) 
Using notations A, B, C, D, we obtain the following relations: 


Sam (8, x) = / ат... ў d, т Чит (hse 25) |. Q0) 


(Messel and Potts 19523) 


O (ни S k\ (t—k\ _ 
Pi (n, теа) = У) == у | es sh neg ә) (1.1) 


t=n+m k=n 


In deriving (1. 1) we made use of the generally known properties of the generating 
function: 


1 oam А 
АС iir en к= ы d 
. . 1 2 


и=из=0 


: Qnam 
Salaa [тат би d 
ш=из=1 
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$ 2. Role of primary particles. Derivation of auxiliary relations 


We make use of the fact shown by Urbanik (1956): 
Sia (5 х) = ете Уз (в, x) 


o CX T^i 0 >0 
AT. B V. (5 x) = M pie (2: 0) 


We insert Si, written in the form above into the equation of the first collision for 
the moments (Janossy 1950) 
2 n—1 m-1 : 
LU , n т 
ее ail Vi „(&®%® = Diu ШІ ІБІНЕ 
Ісі k=1 $ 
1 


XV ink (5.2) Vii bm] + | wi (e’) de’ ie [e 


0 


э ЕБ : 3l (2.1) 


d divide on both sides by pe . We pass to the limit and employ (2. 0). The ratio 


AE - yi жи „m Vanishes, since Vis tends to infinity more slowly than exponentially. 
yi obtain i 
1 
e € — 
a; — 05) = lim w! (ell de | ——————— + м 
Eg ps MM Peru Vi... (6.2) 


And therefore: 
ER 1 


ee). 
0 = lim js (ғ) de’ ле, | 22) 


хо оо 


| 
1 E 2 Ж же. 
04 — Ge = lim CHA OY CES = EK UR TUM евро сте сыш 


x00 Am (e, x) V2 (e, x) 


Using the method given by Stachowiak for the first moments of the electron-photon 
cascade, we obtain, as a result of the fact that relations (2. 2) and e. 3) are identical 


# 
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with Stachowiak’s relations, the following interesting relations: 


: И, (є, x) . о т (=, х) 
pene le ENS = eee — 2; 
lm АЧУ КО ООЛ е0 
1 
E 

eee 
- auc WIESEN (2.5) 

Wa P ТЕС ҮР е Mea URS АГАРЫ 3 

бі — Ga ар ил (el de Е 


(Stachowiak 1956). 


We now consider the conclusions resulting from the equation of the first collision for 
the product densities. This equation has the form 


д 5 > 
(2 + «jen (n, х) 


1 
= M [vo de Е (2.2) dis (ен (2.6) 


k+l=n 0 
p+s=m 


In the notation used, q is the set of variables (hi Nn: Nnti- Пит; п |е is the set of 
variables (mle»--- пъ T iE ж 7 1—2 18 the set of variables (y;|1—e,... 
nill—e; nyill—e,-.. +8), where (т... Nk, т т) are a permutation 
of the variables ur 5), and the variables (Ик... 7,45, Mei +.) ате a permu- 
tation ОШКО. 7„+т). The summation sign with the asterisk denotes the 
summation not only over the indices, but also over all possible combinations 
(71... т), Ope 7), and all posisble combinations (5, 41-2 lad ss (nis Tia: 
This equation is an application to the electron-photon cascade of the equation given 
by Messel (1952) for a nucleon cascade. 

Using the formal analogy between equation (2. 6) and Janossy's equation e we 


can solve this equation formally for n nm m. >> 1 by employing the method given by 
Lopuszanski (1952). 


We obtain!) 
Im sé 
Tu 4 д с 
а Е е Га У wi (2 ж-Әіх 
m 2, 7 ar 7) БЕРІ 55 2 SE (n, х к) 
0 0 p+s=m | 
г As ve 
SC e d д? а) | (2.7) 


! In Eq. (2.7) and after бір 82у is the Kronecker à 


д 


Asymptotic Behaviour of Higher Factorial Moments 309 


The prime beside the summation sign denotes that expressions containing ql,, and 


43-3 should be omitted from the sum. We note that because of (1. 0) and (2. 0) the 


integrals 


Га дь» (т 1) 4857 (о, а)ах; Га д» (ль а) еч д5) (о, x) da 


will exist. 


We shall estimate the ratio 


хус д 
Чат (En) | | 56 CHR (0, ә) 


For this purpose we shall use equation (2. 7) and а relation given by Urbanik (1956). 
On the basis of Urbanik's theorem, the function 


E 9 i і 
pacc = on балдыр (п, x — 5) = v; (n, x — 4) 
is а non-decreasing function of the variable (x—&). Hence 


#(ъх— 8) «iv (m2) 
Thus 


на ОЁ а 
а [ане (t EE = $ fa ІНЕ IUE 
m Е ы prt de 
223 77 баурау (0, а) қалы 5 ате dis" Мя; қ 


р+зз=т 


It is known that 
>, k+p 
a, (258) =0 ta ef — 1 and > 


-r ; 1+; : 
at E AR dro 2 d 1 
ЖЫ Lë 3 :) E 0 if n—® Ж 5 a 


The product 41 bi Gis? = 0 if either of the Om inequalities occurs. If the inequality 
ау + аз > b, + b, occurs, it follows, of course, that at 23 one of the inequalities 
a, > b, ога, > b, holds too; that is why the product CH d will equal zero if 


Kach 1+; 

2 ent 3283 -(—oreo-nu 
i21 

D 
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Since the set of variables (с)... £j, € ... &) is a permutation of (ву... ғ,) and 


the set (күү... аа Ett £j4,) а permutation of (e, 1... € сари) then 


k+p Р [+s " 
> а; + > & = >> ё; 
i= i= 


and the condition for the product equalling zero can be replaced by 


n+m 


02, 8; > 1 (2.9) 


From the work of Konwent and Lopuszaüski (1956) it is known that 


) Dae 9 
lim ШІ 5,5, (n, ТЕ den (0, x) )| 2-10) for у > 0 


EE del 


From the above and from (2.9) and (2.8) it follows that 


l E 49 А п+т 
lim Je, (e, x) 50 Sa, дау (б,ж)| | --0 “for № e > 0 (2.10) 


>00 H 
e i=l 


We shall now estimate the ratio of qf „ to the derivative with the respect to the energy 
of the first factorial moment in a somewhat different way. As known from the theorem 
of Urbanik, the derivative of the first factorial moment with respect to the energy 
tends to zero if x tends to infinity. Therefore we shall be able to choose a sequence 
{x,} tending to infinity such that for every s we have 


9 i д і | 
= ga Sa, jada; (0 ә) < (= 25 Sa, jda; (Ns Xs — D for O& £« x—0 (2.11) 
қ | O> 


Relation (2.11) follows from the fact that the first derivative of the first moment with 
respect to the energy is the product of a decreasing non-positive function argument 
x and an exponential function different. from zero at every point. We have 


ioc д 
ET (=, Xs) (Б 20 58 play (0, all 
2 *;—8 E Ч 
vL IE faj ме B E: Зари D al 2 
jug | 
7 (9 ELT s S Senf nre z 
30 ір E^ » Xs dk,p. P QLs PS , 


k+l=n 
р-Ез-т 
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Thus, from (2.11) 


: 5 > д d 
lim Е (5 3| | ON SEN (0,5) IE > lim у | 085 
DNH 9 А ed Те 
хш” (2) ET An (nad) (35 Shon O29) У x 


k+l=n 
p+s=m 


It is known that 
- k4p 
КӨЗЕ t аш 


>” 1+5 
ЕЕ Meer NS 
and 4, (е) #0 ns <%- о 


From the paper of Konwent and Łopuszański (1956) and from the definition of the 
sequence (x,) 


9 55, us, (р х) 958,8; (1, Xs) 
lim [end b Et == йт Ee «sel So foe жалын i 
E 0) Ss Feo (0 EE CEA 
9 0 | 9 0 
It is readily noted that only for 


k+p 1--5 п+т 


> Dept ep = De; 
ici ігі ізі 


will there exist on plane (e, 7) a domain in which the product 41» (гі = 0 and 
the limit of the ratio of the derivatives with respect to the energy of the first 
moments both tend simultaneously to infinity. Thus, for | 


n+m | 

Da e; lim VIL De NES. drm ES — oo 
0 => 5-СО 0, 

і=1 = San ( Xs) 


Because of the physical meaning of these quantities, this limit cannot depend on the 
choice of the sequence. Therefore 


n+m 


im |4,69 [E 20 БӘ (0, з) |- оо for 0— $ e, (2.13) 


x 


# 
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Dividing both sides of (2.10) by (2.13), we obtain 


д 
| MERC 90 Ee (0, ә) 


Inm (е, а) x lim. А = 0 


lim 


sg - 958 pray (0, 2 ео dos (е, а) 
90 
ndm pts р 
for 0 < 2 e and for 0 > 52 в: hence 
i=l i=l 
(г. n+m Pts 
nm (6, 5). €, X 
lim = 0 for ee (2.14) 
со abe (© - D 23 


$ 3. Conclusions from the equation of the last collision for product densities 
We shall make use of a fact resulting from (1.0), (2.0), and (2.14). On the basis 


of these formulae, the function qf „ can be represented in the form 674% а (в, х). 


The function QU (e, x) tends to infinity less slowly than exponentially. For this 
reason, we may assume 


dic dre 
due Она (ғ, x) 
lim o e же) (3.0) 


х-с Quin (=, х) 


The equation of the last collision for the product densities (Messel and Potts, 
1952b) has the form: 


9 i 
E + n 04 + m м) Qnm (=1... Ens En+1... Ent+m> x) 


i , , ; En 1 
T у G(Qn—2,m4-1 (£1... &n—2; En+1... Entm En + En—15 x) w? Sa arc. i9. 1A) 
En F En—1 


GN En cb En—1 


; F 2 ; e ; 7 et 1 
+ у. А (GE En—1; бөлі ETH EDM Entm—4> хуи кт) “Ks тес (3. 1B) 


GN c, M Ext пт] €n Е Бит 


Ss ee (е1... ыс ев ар (2) | (3.1C) 


cage 


i H , = 
+ 2 ie (61. En U; TS fa i: a) w1 (1 кыш Se er (3.1D) 


c, M 0 


1 
+92 f Qh—1 m1 E E E ёпт, и; X) ий £3 Ga SE) 
CN б 


и 
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Notation: 
NVanables step £,) are а permutation of variables (e, ......... E) 
veriables (6. cum En4m) are a permutation of variables (e, , , ......... CE) 


The summation signs with the indices 1) С, 2) СМ, 3) С.М denote, respectively, 
summation over all possible 1) &,, 2) San 3) е,, Eu , under the condition that if for Dé 
we take є}, then for ғ, 1 we can take only such e for which L+ k. We substitute 
into equation (3.1) the function qf „ in the form of (3.0), after having taken into 


, 


u u п+т 


and the obvious equality 


account the fact that ш! (=) = 0 for и < = and w! ( — 2d =Oforu<e 


Gan (e, x) Бет Оһ» (е, х) 
Чл» (€ х) Оһ, (e^, x) 
we obtain n а (т — 1) а = 


, 


1 , Га ГА , 
2 x lim Lom (Er En—23 Enti Entm Ent E133) 2 = (3.2A) 


CNT фат (Ey Ens Enti- Entms X) [£n + ёп—1] En. En—1 
c дал (1... SCH En EE SEN күз ee x) En 
Ne ee EE EE Lee | 8-2) 
eeh. T Inm (Бү En; Entre: Ent-m3X) [En H Entm] En + En+m 
1 1 | 


1 


TIE Фо (Et Ends 02 Enda: лет) Sn du 9 
TY us [o Qro tecti tern menit | -+ (8.26) 


М x% -+ 00 e Ч Lex End ёп+1" п-т? x) u 
C. € 


1 
ғ“, 


PNE л у ee (3,2D) 
u u 


1 
praia d 8» Qnm (Ey Ел; Ent 1 ит} X) 


1 


MEN S ү ы e ite eat? w? E e (3.25) 


i 2 
Qn.m (ғі-- Ens En+ 1 En+ms х) 


Let us examine Eq. 3.2 for the particular case іп which n — 0. It is readily noted 
that we shall have to discard the following terms: | 
(3.24), (3.2B), (3.2D), (3.2Е). We obtain 
1 


А А 3 , , 2 CL d 
ET һа f nitet us 4 


ТЕР Фб енә) 


314 J. Czerwonko 


For т + 1 the side on the right must be different from zero. From (2.14): 


i 7 , m-1 m 
ЕЕ WIR Bes Bon AA | 
E 1 (us А330) =0 for u+ ej > у =; 
H D 
сс Фот (Ез 6; X) i=1 i=l 
m m 
Since D & = я с; the condition that 
іші i=l 


litio: peal А 


xco 
equal zero is replaced by ғ, < и. Therefore, if we were to have 
1 ГА d , г : 
lun 41т-1 (Ems £r Em—15 X) 
; 
х— оо Jom (£177 Ems X) 


then the integral on the right would vanish and we would have (m — 1) б; = 0. 


mee 


Hence, for m > 1 we have 


SE Gan 3 (Em; Ei ём-1;®) _ 
x00 aon (=1 SE en: x) 


From the paper of Konwent and Łopuszański it is known that for m = 1 this limit is 
equal to zero. Therefore 


оо 


| lm Ят? (Em; бр Em—13 x) 8 irem ли (3.3) 

^ re ДАР, (& эле ens x) ды» 
Let us now pass on to the case in which n + т = 2. We note that by (2.14) the term. 
(3.2C) vanishes. We now write the set of equations of the last collision for n + m = 2, 
taking into account the fact that (3.2C) vanishes with x going to infinity (Eq. (3.3)). 


From the equation for п — 0 it follows, as previously (see (3.3)), that 
: gue (£1, 6; х) 
lim —————— = 
х-со Qi (E15 6 А) 
The remaining equations, after taking into account the above equation, have the form 


do (Es 16 6422) 1 81 | 1 
ш 


04 = lim ———————-— TOC n 3.4 
4 х- Qu СЕ ғә; х) £ 61512182 ( ) 
1 . 
EEN 
+ lim fesses S wt | EJ LT 
хоо.) Qi (E15 693 X) ЕШ 
fen LEE EIER | & | 1 

2 Gy — x, = lim ———— — ——— ët Leen 3.5 
1 a DU gis (eir ға; X) e KEE pes (3.5) 


1 


| | qi (eis и; x) €2 \ du 
Е 
on 579 Too (Е1» Ee? Al иј и 
1 


, 
€2 
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I denote by В; the terms of equation (3.4) not under the integral sign and by В, the 
terms of equation (3.5) not under the integral sign. The following cases are possible: 


Case 1. B, = ат БЕ 0. 

Caso В 0: В, = 204 — Ge 

(n В, ат; В, < 2a, Ge 

Case 4, B, < о; В, < 20 — Ge 
Let us consider them in turn. 


Case 1. This case obviously is not in contradiction with the information in our 
possession, since from the fact the terms under the integral sign are equal to infinity 
at the boundary of the interval of integration we cannot come to any conclusions as 


to the value of the integral. Dividing on both sides the equalities defining Case 1, we 
obtain 


2 peus Е 
Sie А — sl oe (& + =; Hee (£1, ёз; ail 
= Шт 
ey хоо (doi Lë + 6; x)! (91, (2132532) | 
(2 а — 5) w! | ——— 
&1 + 82 
From the paper of Konwent and Lopuszariski (1956) it is known that 
Фо (6 2) 
lim ———— = 
х>со doi (e, х) 
Therefore 
= 4% (ғі, ё; x) | 
lim ————— = 
хоо qiy (£15 £g; X) 
Case 2. From the inequality defining Case 2 and from (3.4) we conclude that 
В 45% (ғі, ғә; x) 
lim = 


———— = we 
х-со Ф (ғ; £3; €) 


Dividing the inequality in the definition of Case 2 on both sides, we conclude that 


lim {970/91} can be equal to infinity; therefore this is not a contradiction. 
2x00 \ 


Case 3. From theinequality contained in the definition of Case 3 and from 
equation (3.5) 

20 du (65 8532) 

lupe е шш 

Spe "Han GE £25 x) 
Dividing on both sides the inequality contained in the definition of this case by the 
equality contained in the definition of this case, we have 

lita dio (GE £55 x) 24 

; x> 41 (£1; £2; x) 

which leads to a contradiction. 


* 
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ee ----- --- 
Case 4. From (3.4) and from the definition of this case 


| que (ғі, 82; X) 
Ша------- 
хоо 411(615 2; 5) 
From (3.5) and froni the definition of this case 
0 44 (E1; €25 %) 
lim ———=— = 
As 4% (GE E23 x) 


which leads to a contradiction. 
Collecting the information from all the possible cases, we obtain 


‚ 911 (E15 6:2) 

lim = = 
X-00 9% (21, £95 x) 

In addition, from Cases 3 and 4 it follows that 


dii (& + 8; х) (201 — gel (ғ + £2) 


х- qao (ғ, E2; х) ipa 21 
& + & 
Employing (3.3) for m = 2 we have 
а (е1, ё; х) : gii (21; £2; x) ж qs (е1, ғә, x) 167; (3.7) 


mo sz a e M mi МЫШ 
x-oo Qi (615 6233) хоо 9%0 (663: X) х-о “20 (E1, Ez, X) 
l 
Doubtlessly, most interesting is the fact that the privileging of electrons for 


n + m = 2 follows from the privileging of the photons for n + m — 1. 
Dividing both sides of equality (3.6) by equality (3.6) taken for ғу and £, such 
that ғ) + ву = &j + & we have 


2 2 
MEL AGO DO О 
lim ST == ————————— 
х-+со 4% (=, 85, х) S | Ei | 
w 
£i + = 


which allows us to generalize relations (3.7). We have 


gho (81, €25%) gi (eiseasx) _ qi, (ei ei; a) 
КЕРИ (САР, eich ae IS ON ҚУЫ f ue DE 
хә Qj (61; &3%) „оо Jao (Er E23 3) — х-оо Чо (Ers &2; X) 


0. (8.7) 


If we investigate on the basis of (3.7) the equation for o m ntm = 3, by 
a method analogous to the case of n +m = 2, we obtain the following relations: 


dis Lex, еҙ, ез; x) E qu (21, 8, ға: Хх) A Gs (ғі; ғә, £3; х) 


dizlbu Gu эсе e E V "m meer ts 
+оо 4192 "5095 63; х х-со 9% (1, £95 53; х) х-со 9 (=, 89% £35 x) 
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qi (ғу, ёз, £3; X) : 81, Eg, 63; 
арта 03 “51 *2» 3 EA Jos (5s 2 А) (3.8) 


х-+со gu (51, 6; 6,1%) асс 930 (61, ёз, 63; X) 


be gio (£1; о» £3; X) 


а +0 
х-со dis (ғі; £5, £35 X) Ge (4 (&, £5, Єз) кар | 


Proceeding similarly as in the case for n + m = 2, we arrive also at the relation: 


M CS £5; X 
2a, = lim HORE UM, а ) ue Le | : + 
хоо ` Ga Lët ғә; 635 X) & + =] et е; 
q^, (21, € + £2; x ; 
E м а) wl | eiie (3.9) 
Cl х-» Qi Lët, 82; 69; х) E2 + ғҙ/ =2- ёз 


1 


We transform this relation. 


204 = lim E ee | b | a ) l SR 


х--со | Q5 (€1, E2; 88; x) & + ё 
34 Е TE (ef, €2, + ёз; ПЕ p1 | E2 ) 1 
x- со gu (е, 81 Sg E25 x). £9 + 23) 5 | £3 


Because of (3.7) we have 


Күн qu (гі, £2 Se £35 д) E? 
x+00 Woo (ғә, 21-Е в»; x) 
from which we obtain 


gi (ез, & + ё; х) 


хоо Why (E1 6; 63; X) 
Therefore 
: Ei, 22 -- Eo, X 7 
204 ren lim Taol ` i w! | Жы. | D (3.9” 
NJ eo Фи (1, E23 ёз; X) 22 F ғҙ/ Ef + з 


1 
From relations (3.8) it follows that the largest among 4*„ for n + m = 3 will 
be qj, since, as may be seen from (3.8), 
NN зв (Ep E» ёз; 5) 
хоо 4% (Ep 82; 83; x) 


for every k= 2; (0 < k < 3). 
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Qi; ou Ee ee ee 
The next largest product densities will be 4% and 41» — their ratio at large depths 

tends to a finite value different from zero. It cannot be determined, by means of the 

methods used for estimation, which of these values is greater. This may depend on 

the specific values of the cross sections. The smallest among the 4,ҙ-ь 15 903: 


§ 4. Conclusions from the relations obtained in § 1, 2, 3 


1 1 
Р 25 пт V] X 
Јат... f d Dm (N> х) (ee (М +m Ze A Ф а) 
e ашы en d. z = 
х- со e 7 ; SH х +00 5 o Q5, n ? x 
Јан... f di. фы (4, х) fan... fangs Soa 
52 г, гі ғ, Qnm (ғ, х) 
1 Pts 1 pts 
n+m Dts ; 
for >) & > >) =, since, on the basis of (2.14), the numerator of the above frac- 
k=1 k=1 


tion tends to zero and the denominator to infinity. For, in the region 
п+т pts i Р 
ЖОО аа 
"em ісі 

which is contained in the domain of integration,.we will have 


lim Zeit _ 2 ER a = ОО 


Уза? бт (Е, x) 


Infpariicular, taking с = $e ина == =... = Ep +s and employing (1.0), 
we obtain 
S (s, x) , 
im gaed 0 for(n+ т) => e (p 4- s) (4.0) 
from which, for & = ғ”, we have 
cue &, X 
m 2%" (50) 0 for (п-т) > (р + 5) (4.1) 


хоо Ba (в, 2) 
In a second particular case: p +s = п +m = М 


Es Ямы (e, x) 


== () for e >" 4,2 

х— оо Shr pp Le, х) ( ) 
Ктот (1.1) 
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Ренье ра дш, LC E senten 
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where R,, is a linear combination of the moments of an order higher than S$ „; 
and therefore, by (4.1), we will have 
ME ies 
lum m ы ( Й )] Е 
х-со Saw (6, X) 


0 


Thus, from the expression Pi (n, m|e, x) and with the help of R,, and Sis we then 
obtain 


By (4.3) and (4.0), we have 


lim Pi (n; т | ғ, x) 
х-со Гі (р, $ =, x) 


= 0 for (n+ m) e > e&' (p+ s) (4.4) 


By a method similar to (4.0) and with the use of the mean value theorem, we shall 
show that if 


lim N^ (59) 0 
soo ам (=, х) 
then also 
lim SN-mm (ғ, 2) Geert 
CO SN (=, x) 
and if / 
lim стт (6 8) _ | At? | 
$ x—+CO N-—s,s (s, x) =Е со 
then we have 
IM SN uM (=, х) x Bet 0 | 


УМЕН СЕ) + > 


By using the relations obtained above, we have from (3.3) 


51-1 (6,4) 1 — бы 


lim — (4.5) 
тека бон (630) біт 
from GA) 
i i i 
EE E (4.6) 
х-+ би (e, x) x00 520 (e, x) х- Szo (в, x) 
from (3.8) 


Sis (=, x) = lim Sio (=, x) -- Jim Sia (=, x) 
x-o біз (s, x) х-о Sz (в, x) х- Sn(ex) 


y 
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о. (4.7) 
x00 521 (=, х) х-+00 530 (=, х) 


im ~ tl 
im WED laaft 


хоо Sh (e, x) 


Employing (4.3), we obtain identical relations for the probabilities P! (n, m|e, x) 


and therefore 


Рі (1, т — Is  l—óy 


ir - = 4,5” 
dun P (0, m| &, x) діл S 
is 210002169) РО Ша) x. P2184) 0 6) 


а О О) аи тері (2, 0| e, x) 


p, DI Ale A wn POIes 1, (1,2,3) 
sues P1216, РОЮ Не: ^ anco Ре 


50 3 | e, x) P((0,3|e x) — 
PIS X P INA) Ue (4.7) 


Ж. Pi(3,0| 2,2) ` BaS] 
mg EE =Е 20 
ES a method similar to (3.6), we obtain 


Jos fii gor ( + е: 4) ШЕ 920 (E, ғә; х) Е ная. E 
q20 (=, £95 x) 


lim = = lim ^ 
х +00 Ska (e, x) Ро 5% (е, х) 


== vmm 


Making use of the mean value theorem іп the integral in the numerator and employing 


(3.6).and (1.0), we obtain 


/ | des qu (E1 + ё; х) 
де) FO 
wl ск cer) 


х-+00 520 (е, х) 


Let us consider the integral in the numerator. Form Abel's Theorem 


1 1 1 1 
К(е) = ] de EK gin (& ++ 2; x) = 2 p d£ ] d£» GO (ғ, + £5; X) 


We employ a change of variables: ғ, + &; = N; dn = dej, and make use of (1.0). 
We have 


1 1 1 1 
В (e) = on q$ (6 -2jas 94 (asd = IE: (2, х) dn 
ғ 26i Р А 


2: 
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Also 
1 
Soa (1, x) +0 
| аера | 
Since 
Sox (7,5) 
lim Зема == (ў for n >> 2e (see 4.0) 


x00 S (=, х) 
then, if we were to have 


Sm (2e, x) E 


lim 
EE 2% (ғ, х) 


we would therefore have 


1 
"EUST ek 
im | SET aso 


x-oco 526 (=, х) 


Непсе 
S Ж (ғ, x) 
lim 


wes SEP 


From the above and from (4.6) 


EN (4.8) 


lim 


PEDE SE (Zera) 


We shall now obtain an analogous relation for n + m = 3. We transform relation 


(3.9). As we know from $ З 


"— 
im Чао (hs "5 А) m+ ni -- 1 (4.9) 
x +оо qag (Ni 172; X) (в) 
т + Ne 


for у + má = т + No 
Taking, in particular, 
4 , 7 ^ 1 
А =: ер таба Nal а ДА а 


and multiplying every term of the sum іп (3.9’) by the left-hand side of the relation 
(4.9), we obtain LEE 
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From the above, as previously, we obtain | 
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Је Ја Јем Jose 5.) 
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then, as previously, we obtain 
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Principle Results 


1) The main factor governing the development of the cascade at large depths 
is the variable formed by the product of the total number of particles and the threshold 
energy. Systems with smaller values of this variable are privileged with respect to 
the systems with larger values of this variable. The primary particle and the qualitative 
composition of the cascade play a considerably smaller rôle these factors are only 
of importance when the two systems have the same value of the “main” variable. 

In particular, for a fixed total number of particles the privileged systems are 
those with the smaller threshold energy; and for a fixed threshold energy, systems 
with a smaller number of particles. 

2) At large depths the factorial moment 52, (ex) tends to the n! m! P' (n,m|e,x). 
This is a simple conclusion from 1). 

3) The róle of the primary particle. Stachowiak's conclusions have been gene- 
ralized for moments of any order. Photons as a primary particle are *stronger^. It 
has not proved possible to say anything about the limits of the ratios 


ТШ: MICE 
im —72————— 
хоо Shy pp (Es А) 
such that 
е ien? 5 (=, x) 
lim ————— = 
хә оо КӘРІС х) 
This is connected with the fact that in the theory two types of equations occur: 
A) those with a fixed primary particle; in this case some different product densities 
for the highest total number of particles enter into the equations 
B) those with a variable primary particle; only one product density for the highest 
total numberof particles enters into these equations. 

4) The role of the qualitative composition of the cascade. А privileged configuration 
of electrons and photons was obtained for n E m = 2 and n + т = 3. One general 
conclusion was drawn for the probability of finding zero electrons and m photons, 
it having been shown that for large depths this probability is considerably smaller 
than the probability of finding m—1 photons and one electron when т >> 1. 


Interpretation of the Results 


We should now interpret a result which at first glance seems to be paradoxical. 
This is the fact that the privileging of electrons for a total number of particles equal 
to two results from the photons being privileged for a total number of particles equal 


ў 


324 J. Czerwonko 


to one. It is however quite clear that if one photon is “stronger“ than the electron, 
then also the “decay products“ of the first photon will be stronger than the decay 
products of the first electron. Thus, the two electrons that result from the “decay” 
of the photon will be stronger than the electron and the photon which are decay 
products of the electron by the radiation of a bremsstrahlung photon. It may be 
foreseen (as confirmed by more accurate calculations), that the “strongest“ confi- 
guration of three particles will be two electrons and one photon, which are the 
result of the radiation of a photon by one of the two electrons of the most probable 
configuration of two particles. 


As to Фе interpretation of relations (4.5) and (4.5’), this is also simple. The 
m photons (m > 1) cannot arise from the configuration of m—1 particles by a simple 
reaction, that is, the radiation of a single bremsstrahlung photon or the formation 
of one electron pair. This fact causes that the probability P* (0, m| ғ, x) is much smaller 
than P'(l, m—1| e, x) for m > 1. On the basis of the examination of the influence 
of the qualitative composition of the cascade for two and three particles, the behaviour 
of the moments of higher order may be predicted. 


It may be noted that the moments of order n +m > 1 may be separated into 
three groups. Each moment oM (e, х) must belong to опе and only one of the three 
groups. This follows from the simple fact, that if we consider a cascade consisting of 
exactly п electrons and m photons, this configuration can arise only by “‘simple 
reactions“ either from primary electron or photon or at all can not arise in апу way. 
By “simple reactions“ we denote.the transitions from occupation numbers (n, m) 
to (n +2, m—1) or from (n, т) to (n, m--1). 

Thus, e. g. 2 electrons and 1 photon arise from the first photon by a sequence of simple 
reactions, since: 

l photon goes over into 2 electrons by pair formation 

2 electrons go over into 2 electrons and 1 photon by the emission of a bremsstrahlung 
photon. = 

Two electrons and one photon cannot arise, e. g. from one electron and one photon 
by a sequence of simple reactions, since this configuration can go over into only two 
configurations of three particles: | 

3 electrons, and 1 electron and 2 photons. 

The division of the factorial moments into groups according to the above criter- 

ion аге; 

l. Those arising form a single photon by a sequence of simple reactions, 

2. Those arising from a single electron by a sequence of simple reactions 

3. Those not being able to arise from any configurations of electrons and photons 
by a sequence of simple reactions — zero electrons and п + т photons for 
n 4 m > 1. 

What belongs the moments of the same group, the ratio of moments for a fixed 
total number of particles and a fixed threshold energy will be tend, for large depths, 
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= 


to a finite value different from zero and dependent on the specific value of the cross 
sections (which we did not consider in this paper). | 

The limit of the ratio of moments from different groups for the same fixed total 
number of particles and threshold energy will tend to zero if in the denominator 
there is a moment of smaller index. 


There arises the question of whether explicit formulae for the configurations 
belonging to the lst and 2nd groups cannot be given. It turns out that this can be done 
owing to the fact that the simple reactions, and hence also a sequence of them, do 
not change the parity of the electrons in the cascade. We thus conclude that to the 
lst group can belong only configurations with an even number of electrons, and to 
the 2nd group, only configurations with an odd number of electrons. Since, as may 
readily be shown, every configuration of an even number of electrons different from 
zero can be obtained from one photon by a sequence of simple reactions, similarly 
as every configuration of an odd number of electrons can be obtained from a single 
electron by a sequence of simple reactions; then to the first group will belong all 
configurations which are even and which have a number of electrons greater than 
zero, and to the second group, all configurations with an odd number of elec- 
trons. 


The properties given above of the division of the configurations of electrons and 
photons may, perhaps, seem doubtful. For purposes of explanation, we add that, 
independently of the possibility of a rigid proof (method of complete induction), 
the properties of the division given are based on a sufficiently clear physical fact, 
namely: if configuration A is *considerably stronger* than configuration B at large 
depths, then the decay products of configuration A will be “considerably stronger“ 
than the decay products of configuration B. This physical fact corresponds mathemati- 
cally to the fact that if configuration А’ arises from configuration A, and configuration 
В” arises from B, where A is “considerably stronger“ than В, then in the equation for 
the “product density“ of configuration A’ will be terms “considerably stronger“ 
than the terms of the equation for the product density of configuration В’. 


We shall now interpret relations (4.8) and (4.10). We assume that we have one 
photon above energye. On the basis of (4.2) we have a photon “infinitely close“ to the 
threshold energy ғ. We now assume that this photon decays into two electrons. Only 
in one case, for electron energies that are equal, does there exist’ the possibility of 
the production of 2 electrons of “infinitely close“ energy &|2. In all other cases we 
shall have only one electron above energy e/2. This was the interpretation of (4.8); 
the interpretation of (4:10) is entirely similar. | 

In conclusion, I would like to thank Doc. Dr. Lopuszariski, under whose direction 
this work was carried out, for suggesting the topic, for frequent consultations, and for 
critical remarks. It should also be mentioned that §2 of this paper is based on the 
generalization and application to the electron-photon cascade equations of a method 
given by Dr. Lopuszanski for nucleon cascades. 
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КРАТКОЕ СОДЕРЖАНИЕ 


Г. Червонко, Acumnmomuueckoe поведение высших факториальных моментов INEK- 
тронно-фотонных ливней на больших глубинах. 


В работе исследуется асимптотическое поведение на больших глубинах по- 
глотитоля высших факториальных моментов электронно-фотонных ливнией 
в аппроксимации „А”. Дифференциальные эффективные сечения — однородны 
и положительны — являются между тем произвольными функциями. При том, 
эффективное сечение на испускание фотонов торможения больше эффективного 
сечения на рождение электронно-позитронных пар. 

Основные результаты: 

1. Главный фактор решающий об поведении ливня на больших глубинах — 
это переменная равна произведению количества частиц на пороговую энергию. 
Система с меньшей „главной” переменной — вероятнее. 

2. Фотоны, как первичные частици — „сильнее” электронов; этой результат 
обобщение на моменты произвольного ранга результатов Стаховяка (1956). 

3. Найдено привилегированные конфигурации системы из двух и трех час- 
тиц. Для двух частиц — это 2 электроны, для трех частиц — два электроны 
и один фотон. 
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In this note a formally exact relation in the field theory is obtained for the elements 
of the S matrix corresponding to two given numbers of particles. All elements with other 
numbers of particles are eliminated. Low’s general equation derived in a previous paper 
(Królikowski, 1957) is a special case, when both distinguished numbers of particles are 
identical. If the two numbers of particles, say №, and N, are different, then the equation 
obtained here connects the transition N, — N, with the transition N, — №, e.g. the 
photoproduction of pions with the Compton effect on nucleons. If the transition amplitude 
for the process N, — М, is known, then that for the process N, — N, satisfies a linear 
equation. 


1. Introduction 


We take into account a system of interacting fields having the total energy H* 
We assume that I) there exist states |A (№) > with definite numbers N of physical 
particles, representing asymptotic forms of scattering eigenstates of total energy of 
the system corresponding to the eigenvalues Е) (№: 


(Ex у — Н) Фм) = 0, 


Фї 1200) + gg; В) (1) 


and II) the operators 
Uy (Е, qr, W) = Py (Н— Е, qo) Put Ux (Р) 
ES (E; (м) UA Py Lu E, (м) Jr ИФ (2) 
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are free from divergence (renormalization), where Py, Ux(W) and V{} are given by 
(5), (10) and (7), or (7’), respectively (Królikowski 1957)? 

We denote by Py the projection operator which acting on the state vector У 
produces its component V^, describing a given number JV of physical particles. The 
projection operator Py = 1 — Py acting on ¥ produces its component Ч, м = 
= V — Wy. N is here (more strictly speaking) a given set of total numbers 
of physical particles corresponding to different fields under consideration. Using 
the orthogonal representation formed by vectors |A(N) >? we may write 


Py= > 1 (№) > <A), “Р = Au Py, (S) 
A (М) N +N 
The eigenstates Ф012, and Ф%7, of H containing outgoing and incoming waves, 
respectively, correspond to the initial conditions У (t) = Wy (to), where tọ > — oo 
and tọ > + co,respectively. 
The time independent equation describing N physical particles (Królikowski 
1957) has the form [W(? (t) = XX? exp (— i W t)]: 


(W — Py AP yp VR 2 = 0, (6) 


if W (to) = Py (t) and to > — >. Here ИХФ is the interaction operator of N particles 
given by the equation 


oo 
Ev ) А Е (+) 
VSP = i [dt Ky (0) ef (NBI + YN) (7) 
0 
2 We note that if 
0 1 
H — H 4 H, $ (3) 
where 
0 
(Exon) — H) 14 (N)> = 0, (4) 
then 


1 
(H — Eau А (N) > = H14 (N) >. 


D а А В d 
We observe that if the conventional renormalization method is applicable H is then identical 


1 
with the (renormalized) energy of free fields and H represents their (renormalized) interaction energy. . 


0 
In the co (3) the operators Py and P LN given by (5) commute with Н. Thus we have 
1 


0 
Д 0 
Е the representation given by vectors |A(N)» were nonorthogonal, then we might define the 
operator Py in (5) using two mutually contragradient systems of vectors (Królikowski 1957). 
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(h — 1), where 
Ree HEIN LN VEEN Pi y HPN: (8) 


We have by virtue of (7) and (6) 


Ee AND (9) 
where 
Dif |! 1 
In (М) = — CIAR Oe о РНР ; 
м (W) ‚| Sak N HP LN 7 Uca mp Nin 


0 


(10) 


For scattering states of N particles equation (6) takes the form (Królikowski 1957) 


(Erm) ER РұНР, SCH US DN = 0 (11) 
and hence 
i 
diy = 1А (М> + ViP(Exan)| А (NI, (12) 


EAN) Ыр р ИЕ 


where VS (Fxg) is given by (2). Here. the solution ФА к of (11) is the com- 
ponent describing N particles of the eigenstate Ф), of Н: 


ФИ м == Рм Oty - (13) 
Note that from the first equation (1) we obtain 


1 


EXE S cw ла Жр HPN OR 14 
Ем) — Pi NHP | N + lE LN М) ) 


(+) 
P | viii == 


if we consider a process for which the system is initially in a state of N particles: 


ШО Ш o) and £y >: оо. 


2. The S Matrix 


Тһе elements of the S matrix for transitions A(N,) > А' (N,) may be represented 
in the form (cf. Lippmann and Schwinger 1950, Gell-Mann and Goldberger 1953, 
Ecksteim 1956, Eckstein and Tanaka 1956) 


€ (АЗ) > = dean — 272 1 (Ем — Baal TYR (ACN), (15) 
where. 
TG. (89) = € AN Бы) 961, >. (16) 


# 
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Formula (16) for scattering amplitudes — owing to (14) — may be rewritten as 


follows 
di GH (4 (N,)) = < A (М;) | Hs N: (H CR Ем.) ie N, D N: (H dd Ew) Р nma di» v 
< ЖЫ, Un, м, (Ns, Ем, Erny) Ф мум, 2 xd 
САМУ, (№, Eray) Banon, » (17) 


where 


Ок, м, (Ni Ем.» W) Tm 


: : Pig HI BONES 
PN, == Fun + (Н Ех) P in, DC DI HD сы NP 
Рм, (H—Eung) Px, i] dt Knie М, Вим f) ets (18) 
0 
(et В 4 GE) 
Vyn, (№, Eina) = Px, (H — Erno) Py, — tf 4 Кум, (№, Ему f) € (PN, HPN, + Fa , 
б г 


апа 
Ky,n, (Ni, Ерм к= Py, (Н — Ехм,) Pin, e PLN, HP LN, Рр HI (20) 


Here by virtue of (6) 


VAN, (Ni, Ern) XN, = Uy, (Ny, Ем» №) DA (21) 
Note that 
Un (Еду, W) = Unn (М, Bee, W), УО = VER (М, Еу), 
Ky (t) = Кум (Ехм)» t). (22) 


From (17) and (12) we obtain 
Tix) AND) = € DON) FA, (М, Ел)! ХУ) > + 


i 
s < A ON 2) | РЫ», (№, Ex). 


И Е 3 (М : 
"per mercem О 


(23) 
Assuming that the vectors Ф» Жыр (X) ате bound states of N particles or 
scattering on bound states) form a complete system in the subspace of N particles and 
taking into account the contragradient system of vectors ^N o ВМ [because 
here the system of eigenstates of (6) affords in general a non-orthogonal basis; the 
"potentials" ИХ? is complex; (Królikowski 1957)], we obtain 
SE (A(N) = < АҚҚ, 2) | Kë UN, Evina) Фм, > 


= UNN LEE, (М, Е) АМ) + 
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СА (Na) | VER, (М, Exon) Фм) (0079? VER, (№, Ey) | (М) 
Ехм)- Erno + ie 


*, 


A"(N,) 


СА (№) | Их, (№, Eun) XN e o. ОХ VER, (М, Едма) | AQ) У. 
+ У : 
bN, 


(Won = Eby — ГЪ). Equation (24) forms a simultaneous system of equations together 
with the equation in which the roles of the lower and upper labels of 6€? and y‘* are 
interchanged and V‘*™ appears instead of VC? (and — i’, instead of +) (Król- 
kowski 1957). 

If the imaginary parts of the operators VS. (Nj) and AERE (N,) are so “small” 
that their influence may be neglected in the numerators of (24), then equation (24) 
takes the following approximative form 


ТҰЗЫ (A (ND) = € A (МЕ, (№, Емо) А (№) + 


AP т м) f 
Жі) Ёхм) — Pany + іє 


2 У «АМ VA, (Ny, Erina) Хм,» <A(NDIV EAN Ем) Aw. GE 
Ем) === Ем, + ily, + ТЕ 


. (25) 


bN, 


In (24) and (25) the first term on the r.h.s. ойтон the Born approximation 


for the transition A(N,) > 2 (М) with the exact „interaction“ operator УС} Li UN Evin, M 


of № particles with N, particles. For М, = №, both equations pass over into the 
corresponding equations (Królikowski 1957). | 
Іп i case №, == №, the equation (25) i EU. be linear with respect to nd „ (0X3) 
PEANG P 25 (A(N,) were known. For ТО Dec (A(N,) we have another NE -linear) 
equation (the one given in paper (Królikowski 1957)), which arises from (25) when 
Be ae №, = N,. This equation does not contain scattering amplitudes other than 
TU d. AN). For example the latter equation may be that for the scattering of 
photons on nucleons (y + № — y + №) and the former may describe photoproduction 
of pions (y + N —> л + №). In this example №, = (Iphoton, 1 nucleon). № = (1 pion, 
1 nucleon) and Tir. (A(N,)) describes photoproduction of pions, whereas 


T Ы Ze (A(N,)) is due to the Compton effect on nucleons. 
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Appendix 


The elements of the S matrix for transitions А(№,) > A’(N,) may be also repre- 
sented by 
(№) (М, > = ema — 27 $9 (Ермә — Eao) Tian Q( (N39) (15) 
(instead of (15)), where 


Т) (М) = € ФАН — Ем) HOW (16) 


(instead of (16)). The amplitudes (16) and (16’) аге, in general, equal only on the 
energy shell. 

The eigenstates dio of H containing incoming waves correspond to the initial 
condition: W(t) = Wy (t) and %-> - oo. Then the time-independent equation 
describing N physical particles has (instead of (6) the form 
[PR (t) = Хх? exp (—iW*t)]: 

(W* — PyHPy — Fei р) = 0, (6) 
where V&?* is the Hermitian conjugate of Vf, that is now defined by the equation 
(instead of (7)) 


eoo (PHR Еу = N) p 
= — if dte Kat (7) 

Ky(t) is given by (8). Here 

УО, ДО = 0% (Fach 
Ux(W) is defined by (10). 
For scattering states of N particles equation (6’) reads 
(Exny — PyHPy — РС») Dw 0, ((11’) 
and hence 


iL 
Pann = AN)» + 


Eis = Pg Pye VR (Ех) ЖАМ)». (19) 
Неге 

ФО = UB NDIN) (13% 
We note that from the first equation (1) it follows that 


1 
Ехму— Рам HP jw — ie ТЕ 


if Y (t) = Fy (to) and г > +o. 


РімбО,- РімНРу Prin 4 (4^) 
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By means of a calculation similar to the previous calculations we obtain for the 
scattering amplitudes (16’) 


Ты) (АО) = <Ф аум, Vv, (№, Ем) ИМ), (170) 
where 


О (PN HP ти); 
ХХ, (Nos Eqs) = Pull — Ej) Py, — i | Шет он NN (No Ему); 
(197 


Кум, (№, Ему t) is given mutatis mutandis, by (20). Here 
NN, (Ns, Е, ом) = ТЕ О.м, (No, E Ny W) KE А (21) 
where Сум, (№, Ем.» W) is defined mutatis mutandis, by (18). Further, we have 
Ta: UU = <Я UE Vi, (№ Бам) A(N) > + 


1 = 
V SN (Na Enn) М, 


ВН Т 2 
all Vn, (Ex) КТС oe an ы: on 


The approximation form of (23’) is 
Tron (А (3) = T (М) V unde Exny)| 5 (№1) + 


2 = ш + ie 
A" (М) 


zd M < Хаң, АМ Excny)| a’ (NA ` GV NON (№, Ex) | XN, )> 
Erna — Eon. tly tte 


(257 


bN, 


The first term on the r. h. s. above represents the Born ра for the 
transition A(N,) > 2' (N,) with the exact ,,interaction” ора VS? м, (№, Ew) 
of №, particles with №, particles. 
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Круликовски, Общие уравнение в теории поля для амплитуды рассеянных частиц. 


В работе получены на основании поля точные соотношения для элемнетов 
матрицы $, отвечающих двоим заданным числам части. Все элементы с другим 
. числом частиц элиминировались. Общее уравнение Лоу выведенное в предыдущей 
работе (Круликовски 1957) есть специальным случаем, в которым два различае- 
мые числа частиц были идентичными. Если эти два числа частиц, предположим 
№ и № с переходом N,—N», н. пр. фотопродукция пионов C эфектом комптона 
Ha нуклеонах. Если амплитуда перехода процесса Мі->Мә известна, тогда ампли- 
туда процесса №№ выполняет линейное уравнение. 


W. Królikowski 
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In this paper the potential barrier height H of internal hindered rotations of methyl 
mercaptan molecules was estimated by means of a method based on the scattering of 
thermal neutrons with a Maxwellian energy distribution (mean neutron energy E = 
= 0.037 eV). The cross section of СН,5Н molecules measured by thermal neutron 
scattering was o = (194.5 + 5) barns. This was determined by the method of relative 
attenuation of the beam. Water was used as a standard liquid (трт у = 91 barns). The 
cross section for rigid CH4SH molecules calculated on the basis of the theory of Sachs 
and Teller was б,;, = 258 barns. 

Assuming a potential barrier height of Hes = (05 а] and taking into account 
the hindered rotation and the vibration of the atoms, the authors calculated the cross 
section of the CH4SH molecules. The result obtained was б), о, = 202 barns. The cross 
section for collisions that are elastic with regard to hindered rotation was also calculated 
under the assumption that H = 1460 cal/mol, this value being obtained by thermo- 
dynamical methods. The discussion of the results leads to the conclusion that a height of 
705 cal/mol found by microwave spectroscopic methods should be ascribed to the potential 

barrier in the methyl mercaptan molecule rather than the value obtained by thermo- 
dynamic methods. 


1. Introduction 


In 1953 and 1954 measurements were made by one of the authors (Janik 1953 
and 1954) of the cross section for diffuse scattering of thermal neutrons by molecules 
of liquid CH,0H. These measurements were carried out in order to investigate the 
influence of hindered rotation in СНОН molecules on neutron scattering and to check 
the theory of this effect as given by Kolos (1954 and 1955). The results obtained were 
within the limits of error, in agreement with this theory, it being assumed that the 
potential barrier height of the hindered rotation іп the СНҘОН molecule had the 
value of 1110 cal/mol (Swalen 1955). 


4 (335) 


356 А. Budzanowski, К. Grotowski, J. A. Janik, W. Kolos 


There exist important arguments, on the basis of which it may be supposed that 
the obtained agreement of experiment with theory was a chance occurrence. It is 
known that methyl alcohol is a strongly associated liquid, which, according to the 
theory of Sachs and Teller (1941), should cause an increase of the cross section for 
neutron scattering as compared to the case with no association. On the other hand, 
however, the discovered vibrations in liquid CH,OH of frequency 250 cm— !, and 
interpreted as vibrations of the hydrogen atom of the OH group perpendicular to 
the direction of the O — H...O hydrogen bond, should produce a reduction of the 
cross section. The agreement obtained between the theoretical and experimental 
values of the cross section might be the result of a chance superposition of the two 
effects. 

This view gave rise to the undertaking described in this paper — the measurement 
of the cross section of the CH4SH molecules for thermal neutron scattering. Because 
of the weakness here of the S — H...S hydrogen bond, the suggestion that the in- 
fluence of the association is not taken into account loses its meaning. 

It should also be noted that the methods commonly used thus far for determining 
the potential barrier height of hindered rotation does not in general lead to consistent 
results. Thus e. g. for mercaptan molecules, Kilb (1955), using the spectroscopic data, 
found the potential barrier H — 705 cal/mol, while thermodynamic methods gave 
the result Н = 1460 cal/mol. 

It might be expected that the results of the cross section measurements for 
neutron scattering in methyl mercaptan would allow the estimation of the potential 
barrier height and eventually decide between the values obtained by the methods 
cited above. 


2. Theoretical Results 


Because of the already mentioned divergence of results for the potential barrier 
height obtained by spectroscopic and thermodynamic methods, the cross section 
of the mercaptan molecules was calculated for both H — 705 cal/mol and H — 1460 
cal/mol. 

For both these values of the potential barrier the coefficients a, of the Fourrier 
expansion for the periodic part of the Mathieu functions describing the hindered 
rotation were calculated! 


Mom (р) Leche > a (eis) ein?” (1) 
n= — 00 
From the formula 
bo (Tis ту) = > а, (о; т) а, (Of ty) | (2) 


+ Рог the notation used in this paper cf. Kolos 1954, 1955. 
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the coefficients bg were then calculated for various collisions taking place with a change 
of quantum rer т (describing the rotation of the entire molecule about the axis 
of hindered rotation), but which are elastic with regard to the hindered rotation. 
Some values of these coefficients are given in the table below. 


Table I 


0€ 1 ТАТ 


Н -- 705 cal/mol 


Н = 1460 cal/mol 


The cross section of mercaptan molecules for collisions that are elastic with 
regard to hindered rotation was calculated from the formula (Kolos 1955)? 


4 Oo (Uo Hy из H3) . DIT AEN RAD a UA 
a ТЕ 2, bà (ti, ту) D ІІ mem (3) 
where 
_ (2n)! (s 4- n)! By’ E; a2/h?)" (d -- m 4- n)! 
E (п!) (n + 1) Y» " (d 4- s+ m+ 2n + 1)! 4) 


The results obtained іп this manner, averaged over the initial states, are shown 
graphically in Fig. 1. 

The vibrations of atoms in the molecule still have to be taken into account. 
Since the calculation of the influence of vibrations on the cross section would be 
very complicated in our case, we have assumed that it is the same percent as in the 
methane molecule (Messiah 1951). s 

In this way, for a neutron energy equal to 0.037 eV, which corresponds to the 
mean value of the energy of a Maxwell distribution at room temperature; the following 
values of the cross section of the mercaptan molecule were obtained 


с = 181 x 10 —?* cm? for Н — 705 cal/mol 
с = 195 x 10 —?* cm? for Ж = 1460 cal/mol 


However, calculations show that for Н = 705 cal/mol the first excited level 
of hindered rotation lies only 150 ст! above the ground state level. As a result the 


? [n the same formula derived in the paper cited, there is a misprint. The results listed there, 
however, were obtained by means of the correct formula whose form is given above. 


228 A. Budzanowski, К. Grotowski, ). А. Janik, W. Kotos 


neglecting of cross section due to collisions that are inelastic as regards hindered 
rotation is not justified. 

By using formula (3) and calculating the coefficients bọ? for the transitions 
connected with the first excitation of the hindered rotation it was found that for 


0.04 0.02 0.03 (ей) 


Fig. 1 Cross section of methyl mercaptan molecule as a function of the initial neutron energy calculated 
by taking into account the effect of rotation and internal hindered rotation. Curve 1 corresponds 
to a potential barrier height Н = 705 cal/mol, while curve 2.to a height Н = 1469 cal/mol. 


Н = 705 cal/mol the cross section corresponding to these collisions is equal to 21 x 
x 10-24, Hence the total cross section of the molecule for this potential barrier height 
1202 Ост? 


3. Experimental method 


The equipment used for measuring the cross section for neutron scattering 
consists of a source of slow neutrons, a scatterer and a detecting arrangement. 

The neutron source was 100 mg of radium mixed with beryllium. The psc 
were obtained by slowing down in paraffin shaped in the form of a ,,neutron gun". 

The scatterer (CH,SH) was placed in a flat-walled brass vessel whose walls were 
‚ thick enough to withstand the vapour pressure at room temperature of the liquid 
under investigation. Because of the fact that the vessel was to be used for investigations 
of other substances with a saturated vapour pressure greater than methyl mercaptan, 
it was made to withstand a pressure of about 100 atmospheres. The brass walls bound- 
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ing the layer of the investigated substance had a thickness of 5 mm for this reason. 
The layer of the substance under investigation had a thickness of 5.3 mm. 

A multi-plate ionization chamber containing a set of 6 electrodes covered with 
a layer of borium of thickness to about 3 mg/cm? and filled with argon to a pressure 
of about 0.5 atmospheres (Budzanowski and Grotowski 1956) was used as a slow- 
-neutron detector. 

The electronic equipment working with the chamber, and designed to detect 
electron pulses arising in the chamber, consisted of a pre-amplifier, a model 100 
amplifier (according to Elmore), a Schmidt discriminator (Elmore), and a scaler. 

Absolute cross section measurements for neutron scattering must be made in 
a suitable ideal geometry in which the scatterer is placed sufficiently far from the 
source and from the detector. However, because of the weak source of neutrons which 
was used in this work, it was decided to work with relatively small distances between 
the source, scatterer and detector; to make up for this, relative measurements were 
made by using water as a standard liquid with a known cross section. In this geometry 
a certain percent of scattered neutrons found their way to the detector. This fact 
caused a deviation from the exponential character of the absorption which was difficult 
to estimate if only absolute measurements were made, 

With the use of relative methods of measurement, however, the matter looks as 
follows: 

The number Г, of neutrons making their way to the detector, with water used as 
the scattering substance, is expressed by the formula 


Ja = Ле e"o + JL — елә") В 


where J, is the intensity of the beam of neutrons falling оп the detector, с, is the 
cross section of the water molecules for neutron scattering, п, is the number of mole- 
cules of water in 1 cm?, x is the thickness of.the scatterer, while J, is a factor which is 
determined from the water measurement and whose value can then be inserted into 
a similar equation written for the scattering liquid under investigation.! Using methyl 
mercaptan as a scattering liquid, we obtain the above formula in the following form 


J, = yen? - Ji (1 — етт") (6) 


Eliminating /, from both equations we obtain 


Т; — Joe w^ ^w * 


1 | Jo. — ПЕЕ дра: (2) 
Zeeche ср ел» 
n ] — egen" "w` * 


EE е. —————э ө —ө————-—— 


1 The factor J, is defined as the percent of scattered neutrons which, because of imperfect geometry, 
enter the detector. , 
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From this formula the cross section is calculated after inserting the experimentally 
measured values x, Jy, І,, and Г,. From the discussion of the errors it follows that 
the errors of the measurement of x and J, influence the value of the error of ош to 
a very small degree, so that the main part of the time of the measurement was devoted 
to the measuring of the values /„ and Iņ. 

The positions of the source, scatterer, and detector in the geometry used for the 
measurements are shown in Fig. 2. 

In order to eliminate any eventual changes with time in the threshold of the 


equipment, the measurements were made by using two identical absorption vessels, 


Fig. 2. — Geometry in which the measurements were made: a) neutron gun, b) vessel with the 
scattering) liguid under investigation, c) 6-plate ionization chamber, d) cadmium shielding againts 
scattered neutrons. 


one of which was filled with water, while the second with mercaptan; the vessels 
were measured alternately over fifteen minute periods. At regular intervals measure- 
ments were made of J, by spilling out the water and using an empty vessel. 

The identicalness of both absorption vessels was also investigated by placing 
them (filled with water) alternately for a number of 15 minute periods in the same 
positions as for the main measurements. It was established that the vessels are not 
completely identical and that the number of counts obtained by using the second опе 
when compared with the number of counts obtained by using the first had to be multi- 
plied by a factor of 0.983. 

In order to calculate п, it was necessary to know the density of mercaptan at the 
temperature during the measurement, that is, at a temperature of 20°С. The measure- 
ment of this density was made on the basis of the values of the density at 0?C as found 
in the tables. This measurement was made as follows: a column of mercaptan of 
a certain height was closed in a thermometric capillary tube. Next the variation of 
the height of the column with temperature was investigated. Because of the consider- 
able distance from the critical temperature, it was possible to neglect the density of 
the vapour. If the value of the density at a temperature of 0°С is known, it is possible 
to calculate from these measurements the value of the density at 20? C. This value 
' was 0.86 g/cm’. 


4. Method of preparing the CH4SH 


Methyl mercaptan was obtained by the method described by Arndt (1921). 
Obtaining CH,SH by this method entails two-steps: First, from a suitable quantity 
of dimethyl sulphate, thiourea, and water there is obtained an intermediary product 
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having the formula (С.Н,М,5)Н,5О,. Then, after boiling this product in a solution 
of NaOH, methyl mercaptan is obtained. The methyl mercaptan is passed through 
a bath of H,SO, and a tube with Cal, 

Details concerning the quantity of substance, the procedure, and so on are 
contained in the above-mentioned paper. 

Condensation of the CH4SH obtained takes place in a spiral glass tube connected 
to the vessel used to collect the CH4SH that has been formed. The tube is immersed 


in finely divided solid CO,. 


5. Results of the Measurements and Discussion 


The cross section obtained in the mesurements for the scattering of thermal 
neutrons by liquid methyl mercaptan molecules is 


см = (194.5 + 5) x 10 -24 emt, 


while the theoretical value of the cross section for the hindered rotation barrier height 
potential Н = 705 cal/mol as given in Section 2 is ø = 202 x 10-24 cm?, 

It should be remembered, however, that the formula (3) was derived for the 
assumption that the energy changes connected with the. excitation of the molecule 
are small in comparison to the energy of the incident neutrons. However, since the 
excitation of hindered rotation is always connected with the excitation of rotation of 
the entire molecule, then in the case of inelastic collisions (as regards hindered rotations) 
the above-mentioned assumption is not satisfied so well. The theoretical result listed 
о = 202 x 10-24 cm? may therefore contain an error. This error would be still greater 
if the barrier height was Н = 1460 cal/mol, since then the first excited level of hindered 
rotation would lie about 210 стт! above the ground state. Without performing the 
calculations, it can be stated, however, that by taking Н = 1460 cal/mol we would 
obtain a cross section for mercaptan greater than in the case of H — 705 cal/mol, 
i е. greater-than 202 x 10-24 em? 

The error which we commit by using formula (3) for inelastic collisions cannot 
really be estimated; it is known however, that it leads to a too great theoretical value 
for the cross section. This follows from the fact that by using the formula we take into 
account also such inelastic (with regard to rotation of the entire niolecule) collisions 
which violate the principle of the conservation of energy. 

It therefore seems that the experimentally obtained value for the cross section 
indicates a height of the potential barrier for hindered rotation H = 705 cal/mol 
rather than the height H = 1460 cal/mol found by thermodynamic. methods. 

A far more certain conclusion regarding the height of this barrier could be 
obtained by using a neutron velocity selector. Then, by measuring the cross section 
as function of the neutron energy, the energies for exciting hindered rotation could : 
be determined directly; an increase in the cross section should then be expected for 


і 
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these energies. On the other hand, knowing these energies, we could calculate the 
hindered rotation potential barrier height without knowing the value of the cross 
section. This would eliminate the effect of the approximations used in the theory of 
scattering. 

The authors would like to thank Professor Henryk Niewodniczanski for his 
interest in this work and for discussions during its course. We thank Dr. Janina 
Janikowa for preparing the methyl mercaptan, and Mrs. W. Szymezykowa for making 
a number of preparatory measurements and other aid. 


КРАТКОЕ СОДЕРЖАНИЕ 


А. Будзановский, К. Гротовский, Я. Яник, В. Колос, Ф. Манявский, 
Г, Ржаны, А. Шкатула, А. Ваниц. Оценка высоты потенциального барьера затор- 
моженной ротации в частице CH4SH при помощи эффекта рассеяния медленных нейтронов. 


В настоящей работе расценена высота потенциального барьера Vo внутренней 
заторможенной ротации частицы метилового меркаптана CH3SH, пользуясь Me- 
тодом рассеяния медленных нейтронов с максвелловским разложением (средняя 
энергия нейтронов Е, =0,037 eV). Измерено эффективное сечение рассеяния Meg- 
ленных нейтронов частицы CH3SH о= (194,55) барнов. Определено оно методом 
относительно ослабления связки. Как эталонная жидкость взята вода Loun —91 
барнов). Эффективное сечение жёсткой частицы CHaSH вычислено на основании 
теории Caxca и Теллера orig=258 барнов. 

Допуская высоту потенциального барьера Vo=705 са/та! и принимая во вни- 
мание заторможенную ротацию, а также вибрацию атомов, вычислено эффек- 
тивное сечение молекулы CH3SH, получая gteor—202 барнов. Учитывая затор- 
моженную ротацию, вычислено тоже, допуская Уо=1460 cal/mal, эффективное 
сечение электрических столкновений. Обсуждение результатов приводит к за- 
ключению, что потенциальному барьеру молекулы метилового меркаптана скорее 
следует приписать высоту VE 750 cal/mal, найденную методом микроволновой 
спектроскопии, чем V =1460 cal/mal, т. e. величину, полученную термодинамиче- 
скими методами. 
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THE APPLICATION OF THE TIGHT BINDING METHOD 
TO THE INVESTIGATION OF ENERGY BANDS 
IN HEXAGONAL CLOSE-PACKED STRUCTURE. I. 


By Marta MrAsEK 
Institute of Theoretical Physics, University of Warsaw, Warsaw 
(Received April 3, 1957) 


Tables for reduction of the matrix components of energy in the tight binding method 
have been prepared for hexagonal close-packed lattices. 


The tight binding or LCAO method originally proposed by Bloch (1928) can 
be applied satisfactorily to the investigation of the symmetry properties of energy 
bands in crystals. Slater and Koster (1954) worked out this method in detail, discussed 
the possible approximations and the validity of the application of this method to real 
cases. | E 4 

А rigorous application of this method is hardly possible on account of the 
enormous amount of numerical work. Taking into account the approximations 
discussed by Slater and Koster, this method can be applied to crystals, in which the 
electrons are tightly bound. The LCAO method can be applied also in the case, when 
the electrons are not tightly bound, but then only as an interpolation method. Тһе 
energy integrals are then not computed analytically but they are used as disposable 
constants which can be found using the energy values calculated by other methods 
at restricted symmetry points in the Brillouin zone. 

The calculations of Slater and Koster have been performed for simple cubic 
structures, namely for a simple cubic lattice; for face- and body-centered cubic lattices 
and for a cubic structure with basis, in particular for the diamond lattice. 

In the present wórk the LCAO method is used to investigate the properties of 
energy bands in hexagonal close-packed structure. The calculations for hexagonal 
close-packed structure are more laborious bacause the unit cell has two non-equi- 


valent atoms and the geometry of this lattice is more complicated, the directions. 


along the axes of rectangular co-ordinate system being non-equivalent. 
There are in literature a few papers on energy bands in hexagonal structure. 
Herring and НШ (1940) have applied to beryllium the method of orthogonalized 
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plane waves. Schiff (1955) investigated the symmetry properties of energy bands in 
titanium by the cellular method. 

The metal to which the LCAO method could be applied is the hexagonal form 
of cobalt; x — Co. Cobalt belongs to the group of transition elements, as Ni and Fe, 
the d states of which can be treated with the tight binding approximation. For in- 
stance Fletcher and Wohlfarth (1951) and Fletcher (1952) investigated with the help 
of this method the energy bands in face-centered Ni. Suffezynski (1956) calculated 
two-center integrals for body-centered Fe. 

Тһе hexagonal close-packed lattice can be described as a simple hexagonal 
lattice with two non-equivalent atoms in the unit cell The primitive translation 
vectors in rectangular co-ordinate system are: A, = (a/2, — үз a[2,0), A, = 
52 (0/23 үз. a[2,0), Аз = (0, 0, c), where a is lattice constant in the horizontal plane. 
The basis vectors are: t, = (0, 0, 0), & = (а/2, үз a/6, c[2), when the origin of the 
co-ordinate system coincide with the position of an atom of the lattice. The simple 
lattice with an atom at (0, 0, 0) is denoted by the index 1, and the lattice with an atom 
at (a/2, V3 a/6, сј?) by the index 2. 

For convenience we use different length units along different axes of the 
rectangular co-ordinate system, namely o = a/2, В = үз а]2, у = c[2 in the x, у, Z 
directions respectively. 

In our calculations we use the model of the hexagonal ideal close-packed lattice, 
i. e., we put the distance between two nearest non-equivalent atoms equal to a; 


I| ==: Then c/a = ж 1.63. This is approximately fulfilled for many 
crystals with hexagonal structure, for instance for 6 — Со we have с/а == 1.62. 


The energy dependence on the wave vector k will be obtained by solving the 
secular equation 


дец (т/п) ua — д E| = 0, 


тп да 

where (mini, are thé matrix components of the Hamiltonian (energy) between cor- 

responding Bloch sums. In the present work these components are computed using 

firstly the general three-center integrals апа subsequently the two-center integrals. 
The general formula for the matrix components of energy is given by 


(тіп), = ехр [ik (t, — GI 5 exp (ikrj) f Pax (pase ы) H y, (r— rj — t,) dr. 


m, n denote here atomic states, 9, are the atomic functions, гу = 1, A, + 1, A, + l; А, 
— translation vectors (l,l, lą are integers), t, — basis vectors. The Hamiltonian 


H =— A + V(r) is an Herniitian operator. V(r) = > (Ue—r) +U (r— ғ — #,)] 
| 1 


is the periodic potential of the whole crystal, thus V(r +r) = V(r). U (r— ғ) is the | 
potential of an isolated ion located at the point гу. The energy integrals (£-integrals) 
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have the form 


Dos Lët T, Ss г тә Ге, (r— t.) He, (r— г — t.) dr. 


The matrix components are calculated for the following states: s, р (3 functions 
of the type x, y, = times /(г), denoted by р;), 4 (5 functions of the type xy, Wize Xe, 
х2—у?, 3z?—r? times /,(г), denoted by d,). Since we have nine such states and 
a unit cell contains two non-equivalent atoms, we construct 18 Bloch sums. The 
general matrix of energy has 18 X 18 = 324 components. 

Using symmetry properties of the Hamiltonian and of atomic functions we can 
perform the reduction of these matrix components. It is convenient at this moment 
to put the origin of the co-ordinate system in the middle point of the line joining two 
non-equivalent atoms in the unit cell. In such a system the basis vectors are f, and 


нао омаи 


2 
The periodic potential takes the form 


"n-X|v(rz-s-5*) &U(r-neg«)] and therefore 


l 
V(r) = V (— г). We also use the fact, that to every vector г; corresponds the 


vector r; = — г. 
We notice the following general relation 


(т/п) Se Г(а/т),. ]* 


Further relations are given in Table I. 


TABLE I 
(s/s) = (8/8) (5/5)ә = (sls)m* 
(spas E —(slpj)sa" (8/р;) 1а = — (slp)n* 
(s/d) = (sl арм” (5/4 е ті (5/4, ТЫ 
Gjlbpui = GjlPjas ` (рур)  GjlPja* 
EI T (p;lPi)22 , (P;IP;)r2 E (р/рда” i$: 
(рада 225 —(p;|dq)22* (РА mE (рады 
(аа ды --(а,/4 )»» (19412 = (аа 
(dild, MI сс (dil адм” к (919, 2 = (d,/d,)ox* gr 


From the total of 324 components only 90 components must now be considered: 45 


components of the type 


(т/п) = д exp (i&rj) T oe (r) H 9, (r— г) dr 


and 45 components of the type 
(тіл = У} exp lik (ry 4-89] | 9) He, (— п — h) dr. 


ў 
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In the above integrals the origin of the co-ordinate system is again located at 
the position of an atom. This system will be used exclusively in the following. 

We make now the nearest-neighbours approximation. Every atom in a hexagonal 
ideal close-packed lattice has 12 nearest neighbours: six neighbours in lattice 1 and 
six in lattice 2. Therefore, in the matrix components Dina) there remain only seven 
E-integrals corresponding to the lattice sites г) (0, 0, 0); (1, — 1, 0), (2, 0, 0), (1, 1, 0), 
(—1, — 1, 0), (—2, 0, 0), (—1, 1, 0) in terms of о, В, у and in (т/п), six E-integrals 
corresponding to м 4- 6: (1, 1/3, 1), (—1, 1/3, 1), (0,— 2/3, 1) (1, 1/3, — 1), 
(0, — 2/3, — 1). So 7 x45 +6 x 45 = 585 E-integrals remain to be calculated. 
Not all of them are independent. The reduction of these integrals must be performed. 
Here we use more particular symmetries of the Hamiltonian and atomic functions. 
We perform the orthogonal transformations which carry the hexagonal close-packed 
lattice into itself. Under these symmetry transformations the Hamiltonian is invariant. 


We take into account the following transformations: 1) reflections in three 
А ; 2 
vertical planes, the normals to which make the angles с л among themselves; 


one of these planes is the plane x = 0, 2) reflection in the horizontal plane z = 0 
From these transformations we can obtain all other necessary transformations. 


Under these transformations the atomic functions transform in some definite 
way. If we take for instance one of the d functions located on the atom ғ), then after 
any one of the above transformations we obtain a linear combination of the d functions 
all located on some other atom, say гу. 


In the first step of the reduction of E-integrals we use only these symmetry trans- 
. formations after which in a given matrix component, say (mo/n,), only E-integrals 
remain with the same index ту and по. By this treatment the non-vanishing matrix 
components can be divided into ten sets. Every set has different relations between 
the E-integrals. Therefore the k dependence of all components in one set is the same. 

First we write down all the copmonents which belong to every set, secondly the 
relations between E-integrals in this set and lastly the А dependence. We use here 


the following abbreviations: & = а k, = ak,/2, n = В, = V3a КЛЫ с == yk = 
= сы? = ak, V (2/3). | 
1. (8/5), (5/7822 — "ha, (а/а), (822 — r*/32? — lu 
Emn (0,0, 0) == 0 
E, (т) — identicalfor six vectors r;corresponding to nearest See RD in lattice 1. 
(т/п) = Emn (0,0, 0) + 2E,,, (2, 0, 0) (2 cos £ cos 7 + cos 2 2) 
2. (5/5) (s|32? — rha (2/2), (322 — 78/322 — т?) о 
Emn (ту + 6) — identical for six vectors (г, + t) corresponding to nearest neigh- 


bours in lattice 2. 
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(minus 2p (»- 5, 1) cos 6 [Е cos Ё cos 5 7 — sin 2 d чв 


+ оен ёна у SE 


З. (5/х)уъ, (Slay) Gy) (2/2 — ур, (9/32 — r2), 
Glocha, (zxz) (мух? — Уу?) за» (xy[z? — т), (ухх) 
Emn (0, 0, 0) = 0 
Eas: (1.0) 

E150) Е, (141, 0) 
РАКЕ (=2 0-0) 

(т/п)уу = 2 [Emn (1, — 1, 0) — Bas, 1, 0)] sin £ sin y + 2i (E, , (1, — 1, 0) + 
+ Epa (1,1, 0)] sin £ cos 7 + Ey, (2, 0, 0) sin 22) 

4. (5/жл»» ($ худа» (х/у)зә» (x — y?) (/32* — he (уху) лә» (21а) ә, Lol — у)», 
(xy[32? — he (уг|хг)у» 


(mjn) = —^E, „ RES біп 6 cos C (sind EL) 


5. ly) GI? — y?) 27) (yx — y?) (32? — r’) (zly) (02 — УЗ — r’) 
E, (0, 0, 0) == 0 | : 
ENTE SE 10 
Ef, 1,0) = E, (—1, 1, 0) 

E, (350,0) =/E, (—2, 0, 0) 
и ОО Сол: ОЕ ПД Оса озу: 
+ 2E, (2, 0, 0) cos 2E + 21 [Ej (1, 1,0) — Enn (1, — 1, 0)] cos & sin у 

6. (shy) (ls — yha (х/х) ә» (lyha Globe Dis — yho (7/32? — r)a 

(lyzi (ауу) лә» (%2/%2) 95 (уг[у#)уз D — Als — Als (х® — 132° — raa 


f 
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1 1 RAT 1 
Ба | 18 ot? 1) = Ей | I | m D EE 1) = Eom | 15, 2 
2 2 
Eon (o SR 3° | = P (о = . * E 1) 


d 2 2 
(т/п): = 2 cos 6 ЕР | ik E : 1) cos 4 cos 37 ER hu -7 l Jesu] 


7 I MAT 2 о 
+ 2i cos С Е cos 6 sin 3) — Ej, (0 та, sin 3 1 


7. (s/z)12 (21322 — т?) 


1 1 D) 
Emn GES = ino (—1 531] = Ena (о, SEH 


1 1 2 
=~ Ean (154-1) =— Fan (—1,5,—1) =— Fe (0,2, -3) 
1 Е =~ ы? 
(т/п) = — 2Е,, ba l} sin б | | 2 cos E sin 3 "—sngm7|— 


i (os ILI 


8. (в/с) лә, (2/2), Disch, (жу Ул» (ахуе (хуГу2)1ә, (х2|х2 — y "Jha, (%2/3z — 1?) 9 


1 
[uidemus eec ги 
1 
= к. ы Pe 


E. (0,— $1) = Ema (0-2 -1}—о 


(т/п), = — ?Е,, | pE 1) sin E sin С (1. 4 + isin 14) 


9. (lich, (y/2)12, (5/22), (уу е» (2/х2 — he, (жу/х2) 19, (у2/х2--у "he, (у2/322—г "he 


1 1 1 
Eos 1, —, 1 = — Dann ar = = => > 
ое 


i-i 
1а gi) 
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(mn) = 2 sinc [= Тон | 1 $ i cos é sin — у + 


2 N P 1 
e Еп (o, қаса ER a? n| > 21віп 5 E= | bu 1) cos é cos E ?] 


2 2 
E НЕ 152 E 
+ Ban (o 3 1) SE n| 
10. (х/х)уу„ СУГУ)» Gool, (2/2), (у21уг) ц, (x? — il — у?) 
аи ‘ae Kee 1,0) = Pu (—1, —1, 0) = Emn (— il; 1; 0) = Бой (b 1, 0) 
Енн (2,0, 0) = Е,,(--2, 0,0) 
(mín) = Emn (0,0, 0) + Aas (1, — 1, 0) cos £ cos n + 2E nn (2,0, 0) cos 2 € 
E-integrals for remaining matrix components vanish for all r, therefore these 
matrix components disappear: 


(s/z) = (s/xz)u = (5/уз) и = (x/z)y = (у/2) = (х[хг)уү 

= (xyz) = (y/xz)u = (у/у2)1 = (z/xy)u = (z[x? — у?) = (2/32? — ғ?) 

= (хујх2) ц = (ху[уғ)ц = (ха)? — у?) = (х2/ z? — hau 

= (уг/х® — у?) 1 = (у2/322 — т?) = 0. 

Now among all the non-vanishing E-integrals only E-integrals for г, = R = (2, 0, 0) 


and г; + t, = 7 = (0, — 2/3, 1) are independent. Therefore, in the second step 
of the reduction, we express all integrals E,,,, (1, —1,0) and Emn (1,1, 0) in terms 


m,n 
of Emn (В) and the integrals E, , (1, 1/3, 1) in terms of Е, (T). Here we use the 
symmetry transformations, which were not taken into account previously. The 
component (ту/по) will now contain the integrals with т and n not necessarily 
equal to ту and nọ. We have only 71 independent L-integrals. The expressions for 


the remaining integrals are given in Table II. 


TABLE II 


1 === 
Еж (l, + 1,0) = > [Esx (В) Уз E.y (Н)] 


JE GE 
E, (£ LO = [+ V3 E (R) — Bay (R) 


P ii = [x УЗ Е, yy: (В) + E, (R)] 
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1 1 
Е, ху (ь4 JE УЗ E, y p SE: SE 
1 LE 1 1 т 
Fonz | 1,——,1 Euro e О 
3 3' 2 
1 = 
Е, уз (1, + 1,0) = ER [= Е, хауз (В) + УЗ Е, ху (R)] 


1 
Е, (1, — 1,0) = a [Exx (В) +3 Е,,(В)) 


1 1 
ps E Al == [Exe (Т) + 3 Eyy (Т) 


1 


Ey y (l, — 1,0) = z [3 E. x (В) + E ту (R)| 
1 
Eyy (. BY i) =4 sae [3 E хх (T) + Е У,у (7)] 


Е,» (1, + 1,0) = + L V3 E, ,(R) — (B). аа СЕ y (В) 


B 
1 
E pes JEE = — V3 [-E,, (T) + E, , (T)] 
1 1 1 — 
Es Е Hi =V3 Ey, ; T. = —-у УЗЕ (Т) 
Е 


ИЕ 2 њиз xoc (В) + Е, ху (В) + ЗЕ, „ау (В). + УЗЕ, xy (R) 


1 1 
pa H y: = = [Es (T) ЗЕ, (Г)] 


1 1 
Pas H 3 : ) 4 [Es (T) + 3E, yz (T)] 

1 1 UE | 
E H т) ма H 2 ЗЕ, (T) + Bue D 

1 = = 
Baas 10) = > [- Erg cy: (0 УЗЕ, yy (В) + УЗЕ, UI + 3E, у (R)] 
1 

Ey у (1.3 Par a) = Ey xy WE JEE USD 4 [-E х,ху (Т) + E Уух? —у* (Т)] 


E т 
Eyr (1, + 1,0) = 2 ГЕ, ae а (В) + V3 Е, Заї--ға (R)] 
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1 E 1 ies 
Äer ap ПЁ es 1] = иЗ Еа ПЕЕ esc V3 E, epe (7) 
3 2 > 
Ej (131,0) = IOS sagt (В) £ V8 E, y (В) УЗ E, (В) — Ey yy (R)] 
1 1 = 
Ey yz baat = 4 [3 Е xz (T) + E, ys (7)| 
poc so = [+03 УЗ Е, ау: (R) + ЗЕ, (R) Ey y (R) F УЗ Ey. (R)] 
TUN EET | 
Ey-y {Lt} = ER [3 £ xx (T) + Ey s, (T)] 
hc! 
- 1 Ж 
By Заа (1, ЗЕ 1,0) = P [+ V3 E dst (R) == Е, 32*—r* (R)] 


1 U 
Е. 1) - "en (к!) о 


12 ives 
Быз (l 10) = >. [Payee (Ю) + V3 EL, 00] 


1 1 Kë d Ax 
Е xz ki JE =. Е PI a) = Let, (7) | "Т 
s | | pg) ORL TT атар Ж” 
ЖЕЛІН Æ 1,0) EZ as (R) SE E, yz (R)| 
P кес A ure o 1 E 


n ate 


Ey (l, — 1,0) = T Bow E (R) + 8: p at Саур 


252 M. Miasek 


Ej уз з—уз (l, — 1,0) = ГЕ [Exi yii: (№) + ЗЕ, xy (Е )] 
i i T) -- 3E T) 
Exi уз хауа 1, a 1 = CH | URS cuu ( ) НЕ XV XY ( ] 
1 1 T 
Е ху ха 1, lj poem (T) + З Еу? уз ( )] 
: ы SH T (T) 
Exy,yz |1 3 1) = Еау |1, n , = ИЗ [Ёш а (T) — Ё, (T)] 
E dn (, 4 19) = УЕ, ммм (R) — Exy i y: (à + ТУЗЕ, a» (В) 
! 1 auc 
Egw-y [51| = УЗ Ia aam (Т) — Es, (T) 
1 A 
Еул (l, 3:10) = > [Esse (R) 4 УЗ Ey. a äen (В) 
H 1 Sta 
Еу |1, pod ? = из ЗЕ. yas |1, үе Epor ES Eys ys, 328 уа (Т) 
Eens (ly 10) = + H Kann (R) — Ез. (В) F — г УЗЕ Ey, ys (R) 


1 1 
Ey, xz h. Ta is E En V3 [E [Е ‘уз. yz Py Exe xz (T)] 


3 
1 
Еу, iy: 1, m. 1 4 1 yz,x* —y* (Т) + ЗЕ, ра (T)] 
Kä и 1 Lu 
Eyes bl] = УЗ зан (1, 31] = - «V3 Esse (Т) 
1 Ж 
Ел уз зза (l, + 1,0) = = LE äere (R) E ИЗ Ej äer un LÉI 


И 


We could use for R and T instead of (2, 0, 0) and (0, — 2/3, 1) some other г, and 
г, + t, given above. We choose these particular В and T because the direction cosines 
(1, m, n) of these vectors are given by the most simple numbers possible. The matrix 
components expressed in terms of the three-center E-integrals in the nearest-neigh- 
bours approximation are summarized in Table III. 


(sls) 


(5/5) 2 


s[x) 


(5/ж)» 


(s/y) 


Lech 


(5/2)у» 


Leo, 


(s[xy)ia 


(slyz)i2 


(5/22), 
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TABLE III 
Е, (0) + 2E, (К) (2 соз cosy + cos 2£) 


А 1 2 1 2 
2Е (Т) созё [Е ов 7 + cos =) + i ( ion 7] — sin SÉ 


— 2V3 E, (R) sin E sin a + 2 i E, (R) (sin Е cos + sin 2 £) 


\ 


е 1 1 
2/3 E, у (T) sin 6 cos б [s 4 1 — icos Сый 


Е,,(0) — 2Е,, (В) (cos соз!) — cos 2 5) + 2/3 i E,» (В) cos 6 sin n 


1 2 1 2 
— 2E, , (T) cos 6 Ie E cos E: 7 — cos i +i c Ё sin Е: 7] + sin iJ 


1 2 1 2 
— 2E, , (T) sin C IEEE —i ( Е iie. jJ] 


АРИУ Е, x? — уз (В) sin & sin n + 2 i E, „„(К) (sin & cos 7] + sin 2 6) 
— 1 1 
2/3 E, —y:(T) sin Ё cos Ё [o тү 7- us n) 


2E 7) si S sj : mes с & | 2 1 
ы ---- — — $ — -- со — 
әуе (2) sin 6 aues die 7 i | cos & co SS E 7 | 


ЕЗ 1 1 
2 У3 E, у (Т) sin ё sin 6 (s re 7 + i sin Ei n) 


(slx? — y?) Е, хз — уз (0) — 2Е, уз — у: (В) (cos ё cos 7] — cos 2 6) + 273 1 Е, xy (В) cos 2 sin у 


1 2 
(5/42 — y?) mU m -у(7) cos 6 ГЕ € cos sad — cos E n) + 


D D 1 . 2 
БЕЗІ и 


45/322 — r?) Е 322° (0) + 2E, зг: _„ (К) (2 cos 6 cos n + cos 2 2) 


1 2 
(5322 — т?) 2 E, зз, (T) cos б |: cos & cos E 7 + cos = 7 + $ 


(“хп 


(x/x)12 


1 2 
; [2 29 M nep I qud; 
E сов é sin 5 1 m 2] 


Ey x (0) + ТЕ, (В) + З Е, y (R)] cos © cos n + 2 E, (В) cos 2 E 


1 2 
сов б кыт) + 3 E, „ (T)] cos 6 cos EES n --2Е, ДТ) cos de d gp 


1 202 
+ i соѕ б (т kä. Ep (TD) qos sin. n= 2 Exx (T) Ss j| 


х 
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(y УЗ [Ey (R) — Ez, (R)] sin £sinz — 2 i Ey, y (В) (2 sin £ cos ņ — sin 2 6) 


it 
«уу V3 [E, хх (T) — Ey, al )] sin £ cos С c = 1] — i cos тҮ 7 


-— 1 1 
(2/2), 2V3 Ey x (Т) sin £ sin € c 7] +i sin E n) 


(xlxy)u 


E, «,(0) + ТЕ, (НЕ) + 3 Ey a уа (R)] соз £ cos n + 2 E, yy (В) cos 2 E + 


V3 i [Exx (В) + Ey жу (R)] cos £ sin ; 


1 2 
(х/ху) SH se (7) + 3 E, уа (Т )] cos & cos = EE (T) cos El d + 


1 2 
ties с ЦЕ T ates Es cl = ST) sin n] 


шей 1 1 
(xlyz)g = (7127 V3 [Ех (T) — Е’ у; (T)] sin Eent (2 nti sin — d 


(2/22), 


d Ge AT 


(х/х2-у?) 
= Dachs 


(x/3z*—1?),, 
(213z? — г), 
(iy) 


(lys 


Dich 


1 2 
sin 6 E Ey xz (T) +3 Ey, уг ( Т) сов £ sin e 2 E "T -- 


1 2 
+ isin SIM (T) + 3 Ey yz (T)] cos & cos py 28 9) Ser d 


V3 [Ey y? (В) — Ex „у (R)] sin 6 sin — i { [Е x,x2—y? (В) — 


— ЗЕ, ху (R)] sin £ cos n — 2 E, a um (R) sin 2 2) 
il 1 
V3 [Ex xy (T) — Ey.» (T)] sin £ cos х (s С n) 
— 2/3 Е,, 311 — r: (R) sin 6 sin 7 + 2 i Es „а (В) (sin £ созу + sin 2 2) 
1 1 
2V3 Ey 32: p: (T) sin £ cos da SS i cos Е. n) 


Ey (0) + [Ey, y (R) + 3 E, (R)] cos £ cos n + 2 Ey, y (B) cos 2 E 


1 2 
cos С Ж (T) +3 Exx (T)] cos & си + 2 Е, , (T) cos s | 2 
1 E 
+ i cos C 1 [Ey, , (T) + 3 Ez, (T)] cos £ sin — utet Е,» (Т) sin ХО? 


E 9 
2 Е, , (T) sinc Я = 


Ж ER MET 2. 
0 cos ¢ cos ы age 


(Doha 


(у/уг)\» 
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са 
Cu 
On 


УЗ [Ey um (В) — Е, лу (R)] sin £ sinn — 


— i ШЕ. (R) —3 Жоор чур (R)] sin é cos 7 — 2 ET (В) sin 2 =} 
2 ў m an 
йй б FAT A s ЗЕ, xz (T)] cos sin = N + 2E, у: (Г) sin z nj + 


1 2 
+ тт ё | Гры» Е ЗЕ (1) |6 Е о --2Е,,,(7) cos E | 


Mus 


(yl? — y?) BOO [ES SR) с ЗЕ, ху (R)] cos £ cos + 2Е, x1 уа (В) cos 26 — 


Van [Ex «y: (В) + Е, xy (R)] соз £ sin n 


1 2 
(yix? — y?) cos Č | [ Ej (T) + ЗЕ, (Т')] соз.ё cos "RÀ + 2E, à ys (T) cos Sr d ES 


А T wes 
+ і соѕ С E (T) + 3E, xy (T)] cos £ sin "i 7] — 2E, us (T) sin eg | 


(y[3z? — r*),, Ey 32:1: (О) — 2Ey 34: ,: (В) (cos 2 cos n — cos2£) + 2V3i Ex 322 уз (В) cos 6 sing 


1 2 | 1 2 
(y/3z? — т?) = 2Ey 322-12 (7) cos С IE é coe 205 - = a (= é oed] + nan) | 
(2/2 Е, (0) З 2822 (В) (2 cos € cos 7 + cos 2£) 
1 2 1 272 
(212), 2E, z (T) cos C г cos é cos Е. 7 + cos ES d +i ( cos é sin zs 7] — sin E d 
= 1 sc 1 
(alxy); 2/3 E, y (Т) sin £ зіп С (e E n + i sin Ss n) 
(zlyz)ı Е, yz (0) — 2Е, уг (R) (cos £ coss] — cos 28) + 2V3 i E, xz (R) cos E sin y 
1 2 А кы E? 

(2/у2719 —2Е, yz (T) cas ¢ | | cos £ cos — n — cos — n | + і | cos 6 зїп — n + sin — 9 

Б E 3 3 3 
(2122) — 2 ЗЕ: уз (R) sin E sin n + 2 i E, xz (В) X (sin 6 cos 7 + sin 24) 

Из m. s geg 

(2]х2);5 2V3 E, уг (7) sin & cos £ | sin л! — i cos = 7 


П 2 1 2 
(zl? — y?) 2E, am (T) sin ГЕ é sin 22 n+ эшең n) = (> & cos nd — cos En d 


1 2 
{28° — r?) - e 2E 3z: yi (2 ein d (6 сов 6 sin ue sin is) = 


КЕ t 1 472 
LL cos Cos co Te E CUOI. d] 


Е 
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adi NENNEN "== ee 


(xy|xy)i E а (ӨЛЕР Елу, (B) + ЗЕ а уа уі (R)] х cos £ cos n --2Е,,,, (R) cos 22 
Е L ‚ 
(ху/ху)р cos б MES y (T) + 3 Easy (T)] cos 2 cos Es n + 
2 
IE с: (T) cos = ne + 
1 
+ icos С YES xy (1) + 3 Ep yy : (T)] cos б sin Es Я) = 


2 
—2 Е уж (T) sin 2a 


(xyly2)12 


1 BENE 
а V3 [Ey жк (ТА Eye y: (Т)] sin E sin ¢ (> a 7 + isin = n) 


1 
Losch sin б | E [Е ху,х2 (7) +3E ‘yz,x*—y? (T)] соз Be Gp ap 
2 
+ 2 Ё x2 (T) m d $ 
1 
+ isin€ (6 de (T) + 3 Ey, y: (T)] cos 2 cos ТТ n+ 


2 
T 2B ye (ТГ) cos zs d 


(ху/х® — y?) УЗ [Ey ys (В) — E, уху (В) sin P sin 7) — 
— 2 i E, y, y: (В) (2 sin £ cosy — sin 2 ©) 
(хуја? — у) ИЗ [E,y y (T) — Ел ys y (T)] sin £ cos x 
KH к 1 
x cuota ge ет) 
(ху/34® — г?) 2-2 V3 3 Ey ys, zen —— (В) sing зїп] + 2i E xy,3z'—ri (В) (sin £ cos 7 + sin 22) 
Seer Vs 1 d 
(xy/3z? — т?) 2V3 Ex: 2342: 5: (Т) sin & cos C (sin SC 7] — i cos zy n 
Gzlyz)u Ey, yz (0) + [Ez yz (R) + ЗЕ, ха (R)] cos & cos n + 2E ye, yz (В) cos 2 $ 
1 
(уг/уг)уә cos С les (7) + ЗЕ, uz (T)] cos $ cos 3 7-6 


2 
+ 2E уз, yz (Т) cos 3 7) | + 


EPS тысы» — mee 


yz|xz)y 


(уг[хг)» 


(yz/x? — у?) 


(2/32? — rie 


ELS NS dem Pi pigs gus à. = 


ET 
rae “ыз 
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у 1 
+ i cos idis... (T) + 3E es (Т)|сов & sin Am 


ee 
ТҰ 2Е, yz (7) sin 4 7 


е 


УЗ [Ey (В) — Е, (К) | sin Es sin n — 


— 2i Ir (R) (2 sin £ cos  — sin 2 5 


1 
Уз [Е М tos EIER 


Sr 1 
sin 6 | == а (Ту ЗЕ y, y xz (Т)] cos £ sin 574 + 


е 


2y 
+ 2Е, ха y? (T) sin з | Æ 
i sin б lu әкесі D) s S, ух Т)]соз Ё cos ED i 


2 А 
qs ККЕ (Т) сов Si 


: | Й d ` 
2Е, 2.321 (T) sin б Б é sin EN 7 + sin SS 7} | E 


zt. 1 іл Гл» 
mm pmi] 4 


prc 2 3E СРД cos 26 


Sorts “5 Wi Vi 37, E j gd 


Lë 
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1 
(х2 — y?[x? — y?) соѕ б и (T) + 3 E, ху (Т)] cos & cos 2 quae 
2 
+ 2 TURPE ау; (T) cos 3 7) == 
u . 1 
+ i cos C (Еа. Lus ys (T) + 3 Exy ,,(Т)) cos 6 die 9 — 


2 
— 2 E уз ys t yt (T) sin A j| 
(22 — 9/37 — rh. Вл ая (0) — 2Ез rg (R) (cos & cos 77 — cos 2 2) + 


4 2V3i Ey 32252 (E) cos € sin 7 


| 1 2 
(х2 — у2[322 — г), —2 E уз, Зза ра (T) cos e| (= & cos г. 7] — cos m n) + 


1 2 
Je n Lego Zu tn л) 


(322 — 77/32? — r?°)ı Esos sas „в (0) + 2 Ез, зз: (Н) (2 сов 6 cos 7] + cos 3 2) 


ү 1 2 
(322 — г2[322 — 1), 2 Ege 523; (Т) cos t | ( cos & cos в” + cos 2) y 


bod Той 
di oW анасы sure tur elt] 


The E-integrals for r, =0 can be written in the following form: 
Kall = f Pn (® Ho, (r— г) dr 
= f on (9-4 +U (P) +U (r— гу] e, (е) dr + 
tX Tei 0(е-ғ)ө,(ғ-ғ) dr+ 
1-1,0 


пр, ] Pm (ғ) U (r— rj — 5) Ф, (г--г) ағ, 


in which the second and third term contain only three-center integrals. We now make 
the approximation of neglecting all three-center integrals. The E-integrals for г = 0 
in the two-center approximation can be written 


Е„„ (0) = f e (0 E A U() + DUG г) LU (r— ri — &) e, (r) dr, 


where ғ; and г; + t, are here the position vectors of the 12 nearest neighbours. 
Slater and Koster have given a table for the reduction of E-integrals to two-center 
integrals calculated for functions with o-, л-, Ó-symmetry. In the two-center 
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approximation only 16 independent integrals remain: sy, ру, do, d,, do, (ра), (sso), 
(spo), (sdo),, (ppo), (ppm), (pdo), (pdz),, (ddo),, (ddx),, (ddó),. Using Table I 
of Slater and Koster we must know only the directions cosines of the vectors R and T. 
We have for R: 1=1, m=0, n = 0 and for 7: 1=0, т = — V3J3, p= V2/3 


The E-integrals in terms of the two-center integrals are given in Table IV. 


TABLE IV 


ES (О) = Ez (0) = Е, за, (0) = Еу (0) = Ej за, (0) = Қа yas (0) = 0 
E, y (№) = Е xy (R) = E, y (R) — E, ұу (R) — E уха y: (R) — EROS (R) = Ег уз (R) 

= Еул у: (R) = Ey (В) = Ej, ,, (В) = 0 
E, ; (В) = so 

Ех (0) = Ey, (0) = Ez, (0) = p, 

Exy,xy (0) = Ел. зе. уз (0) = do 
Ex, xz (0) = Е,, уг (0) = d, 

E32, 572,322) (0) = d; 

E, xy (0) = Ey (0) = (pd), 

Ej, (В) = E, (T) = (sso) 


2 ít 
E, x (R) — —V3 E, y (T) = үз Е, = (Т) = (spo), 


E Ф) = VBE. 


Sy 


3z?—r? (T) Cd Ey 5,32 —r? (R) = — V3 E 


$,x5—y 


1 jar 1 
M ssm- len 


Ex (R) E (pp о), | 
Вых (T) = Е,у (R) ae Ezz (R) = (ppm), 


$,x1—y* (R) = 


1 2 
Ey, (1) => (рро), + = (ррл), 
2 1 
Е,.(Т) = T (ppo) + T. (ррл), 
Т І V2 = V2 ( л) 
Еу а ) = – 3 (рро); + З PPT); 


— 3 
E, ха (В) = — V3 Е, ху (Т) - Vi Е, ха (Т) = Еу,» (R) — (рал), 


F 
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1 /- 1 
E, ape (R) = — = УЗ Еау (R) = — > (pde) 
"UH 1 2 
Bue (T) = = V2 (pdo), + гү: (рал), 
1 2 
Е, қа. уа (7) = c Meus AF 573 (ріл), 
1 /- 2 
Буза (Т) = — = 03 (pdo), + — (рал), 
2 | 1 
Буе (T) = — 5 (pdo), + 2-03 (ріл), 


es а 2 
Е, а-у (T) = — = 02 (pdoh + Из (рал), 


1 F TE 
E, ззат (T) = ә V5 (рас), + Zr V2 (pd 7), 


> 


Eves (R) = (ddà), i MS 


it 2 { = 
Eo c ола ед a те x 
2 1 Es = 
Livy өг) D = д + = (dd6), еч M 
| О т] E E CECR i MD v 
Жз» (Т) = 5 Cdo + dda) сад} e гй мый = АШКЫ 


"rcp Ss Pe = 
Қас ані n eie o 
SN SE EE Wi v 


АР EDO Ga > 
T DH " 
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Ba äer (Т) = — Vs (ddo), + EN (аал), + E (ddó), 


Ey yas та (В) = — — St (ddo), = — sn (446), 


ie 2 Scis 
Ej ace (Т) = — 1 V3 (ddo), 4- E V3 (айл), — e УЗ (446), 


The table given by Slater and Koster simplifies also considerably the calculations 
of the overlap integrals: 


f фи (r—r;) P(r — г) dr. 


We summarize finally in Tably V the matrix components of energy in the two-center 


approximation. 
TABLE V 
(sls) So + 2(55 о), (2 cos & cos 17 + cos 2 2) 
1 D 
(5/5) 12 2(55 с), cos d ( cos é cos E 7 + cos ES ) + 
JE NL ‚ 2 
+i SE 
(slx)g 2 i (spo), (sin 6 cos n + sin 2 4) 
1 QV T 
(5/^)12 — 2 (spo), sin £ cos 6 (s 50 т ) 
(sy) 2V3 i (spo), cos 6 sin 7) 
2 1 2 
(5/У)12 ue (spo), cos 6 IC Ё cos A nN — cos Sg n) + 


1 2 
+i СЕС 


2 1 7442 
(5/2) - 2 — (spo), sin 6 2 cos Ё sin — 7 —sin — 7} — 
3 2829) 3 
als 1 2 
—i DRE OS а c 


(s/xy) 1 —3 (sdo); sin E sin 7 


1 1 
(5/ху)ю — (540), sin $ cos 6 (o zt n — i cos 3 d 


к 
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(slyzha 


(5/х2) |» 
в у) 


(з/х® — у?)» 


(5/32? — ru 


(s/3z? — г?) 


(%/*)11 


(xx) 


Loic 


(x/y)12 


(x[z)1 


(х/ху)уу 


Locke 


pl 1 2 
Уз (sdo), эт 6 (= F 7j + sin - n) — 
1 2 
= —7— LI 
i | cos 6 cos 3 7] — cos 3 ] 


1 meet 
— 2 V2 (sdo), sin & sin 6 CE Лат QI y E 


— V3 (sda), (cos £ cos rj — cos 2 2) 


l ,— 1 2 
zt V3 (sdo); cos С cos & cos E n — cos СУ n] + 


1 2 
. Pn Lo М су, Ж 
I 3 7] t sin 3 n) | 


— (540), (2 cos ё соз 7 + сов 2 6) 


1 2 
(sdo), cos 6 [Е cos Ё cos En 7 + cos = n) + 


D 2 E si 1 . 2 
+i | (2 сов аш Rea a) 


Po + [(ppo), + 3 (ррл),| cos 6 cos n + 2 (pp a); cos 2 6 


П 2 
cost ((рро), + 3 (ррл),| cos & cos zd + 2 (pp 7), cos i d ae 


1 2 
+ і соѕ 6 | [(ppo), + 3 (ррл),| cos & sin SE: 2 (pp л); sin тұ d 
--Ү3 [(ppo), — (ррл) | sin £ sin o 


1 — 1 1 
=. V3 [(рро), — (ррл),| sin £ cost (s = n — i cos m | 


2 1 1 
59 V5 [(рро), — (ppz)] sin 6 sin С (» ex + i sin E 


2 
(pd), + i E (рав), + уЗ (nt, | cos & sin 7] 


1 | Е 
сов 6 IS (раб), + 5 3 (ез), | со5 2 сов Zn = 


+ 


2 V3 (рал) 2 
3 рил) E, 


i; 1 
+ і сов С Б (pdo), -- 5 3 i, cos Ésin — n zr 


20 ee 2 
Р 53 (рал), o Zi 


(xlyz)ia = (ylxz);s 
= (z|xy)a 


(x[ xz) 


Gi — y?) 
(ха у RÉI 
= (ylxyha 

(x/32? — "Уі 
(x/3z? — г?),, 


NEID 


% Т 


(y/2)12 


(vlyzha 
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D 22 1 
- ГЕ (раб), — EQ? (рал), | sin 6 sin Ё [s E: 7) + i sin i 
: * E 1 
sin б | Е D (pdo), -- 2 y 5e» Е © sin TE n+ 
2 A d. yee 
= ae умы а» 
| E 1 
+ ising {|v (pdo), ДР 2 Уз ca, | cos 6 cos E n+ 
2 2 d. 2 
SR = (pdz), cos xi 


T Ce e = 
Е Е УЗ (рад), — sean | sin £ cos n + УЗ (pdo) өзі) 


Gi 


ВИ = 1 1 
г = EZ (рас), + 5 EM sin 6 cos б [ss = 7 — deem n) 


— i (pdo); (sin € cos 7 + sin 2 8) 


ec 1 1 

— | оз, —— 125 EN sin & cos 6 [is — N — і cos — n) 
3 3 3 

Ро + [3 (рро), + (ррл),| cos & cos n + 2 (pp), cos 2 € 


1 1 
m cos Č | [(рро)у + 11 (ррл)у| cos & cos’ Г) n+ 
2 
+ 2 [(ppo), + 2 (ррл),| cos i л че 
pi | | an 
+ СЕ i cos C. l [(ppo), + 11 (ррл), | cos 6 sin rs 7- 
2 | Е 
— 2 [(рро), + 2 (ррл),] sin = d 
2 1 noo 
EG V2 [(ppo), — (ррл),] sin ¢ | e E n) - 


ver 1 
— L| cos т ы чел 7 


і ЦЕ үз (рас), — (рал), IE Ecos 7 + 2 (ріл), sin 2 д 


= 1 
= Е |" (pdo), -- 10 хее | сос э] + 


D 
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2 2 
SCH ЕН (pdo), + Уз can, | sin ij E) 
1 = 2 1 
E ту isin C 1 (pdo) + o/s cram | cos é cos m + 


S 2 2 
Le Е (pds), + = iam, ei 


3 = 
(yl — у?) (pd), — i 5 (pdo), + УЗ рал), | cos È sin 7] 
(y[a3 — y?) LE m - Тұ ал) АЖА ақы. п + 
7 yu 3 3 MCA MS ролл 3 
2: (54 2 УЗ (pd : 
SE 3 (рас), + 3 (рал), poe ra net 
n T ТУЗ (pd ue 
ck aoe o us 0); — з (рал), eee ys 


Жаы 2 V3 (pd tm 
== ray 0); T z (p л), SE 7] 


(у/322 — r? —V3i (pdo), cos 6 sin 2] 
1 - 1 2 
(у[322 — г?) S |" (pd о), — 4 (ра d cos C | (cose cos E 7] — cos ES n) + 
А Esi 1 x 2 
+1 | cos ry ud эша? 
(2/2 Ро + 2 (pp л); (2 сов 2 cos 7 + cos 2 5) 
2 1 2 
(2/2) Y Ё (рр о)» + ech cos 6 |: cos 6 cos ER 7 + cos E n) + 
Не 
i | 2 cos ani 
(zlyz) 2V3i (рал), cos Ё sin N 
2 L= 1 
(zlyz)a — |: (pdo) — — V3 (рал), | cos С Ic Ё cos — 2] — cos EZ n) + 
3 3 3 3 
D . 1 D 2 
Sg GE 
(zlxz)g 2 i (pdr); (sin ё cos у + sin 2 £) 


2 = 
(zza ~ pis [2 V3 (рас), — (pdz),] sin £ cos 6 (s $ 7] — i cos = d 
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1 | - m 
(z/x? — y?) — — Е (pdo), — 2 V5 (рал), | sin С ( сов £ sin B n+ ee 7 - 
3 3 \ 3 3 
Lien 2 оов — ч 
1 COS С cos E 7) 5008 ту 7 
2 b 
(21322 — т?) zs ПЕ (рас): + тү V2 EM эт С [Е соз Ё ѕіп I + sin n) = 
fe Е 1 2 
—$ cos С cos eg A) Ae GE Yi 7 


9 3 
(хуу) dy + IE (ddo), + (ddr), + ES «un, | cos Ё cos 77 + 2 (дал), cos 2& 


1 
(xy/xy)12 cos 6 = [(ddo), + 4 (дал), + 11 (ddé),] cos £ cos = п + 
2 2 
+ т [(ddz), + 2 (ddó)] cos Ss j| + 
+ і соѕ Ё t [(ddo), + 4 (дал), + 11 (446),| cos 4 sin = n— 
2 2 
EE. [(94л), + 2 (dd9),] sin = d 


2 1 4 11 


1 WES 
x neri den eH 
€ 1 1 
(xy[xz)19 — V2 sint г (ао), — (dd5),] cos & sin a n+ 
2 2 
am TA [(ddz), — (dd6),] sin T d + 
T 1 1 
+V2 isint | [(ddo), — (dd9),] cos $ cos ed 
А аад x 
== 3 [(ddz), cM IN cos 3 Ui 
Gayla? — 9%. VB В (ddo) — 4 (dda), + (dd) sin £ sin 
(ху]х® — у?) Е EEN E (ddo) Sei MS SECH (ddó), | sin £ cos ¢ x 
а 3 4 1 3 4 5 1 


EE ) 1 
х an ат deeem d 


13 
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3i 
(ху[322 — г?)ц E [(ddo), - (dd9),] sin £ sin 2) 


1 4 5 к А 
(ху|322 — г?) — E (440), — E (ddn); + yz «un, | sin 6 cos € X 
D 1 D 1 
X pz piacere 
(у=/у=) 1 d, + [3 (дал), + (dd9),] cos £ cos + 2 (ddó), cos 2 6 
1 19 1 1 
(у21у2),ә а G 2 (ddo), + с. (ddr), + SC (ddó), | сов £ cos a n + 


1 2 2 
dr E (ddo), + E (ddz), + EF aun, | cos = 7 | -F 
1 19 11 RE 
+ = i cos G 4| 2 (ddo), + Es (ddr), + Ss (ddo); | X cos E sin SS n — 


1 2 2 
12 | (ddo), + —- (dda), + — aan, | sin j| 


(У2122)11 = (ал), - (ddó),] sin £ sin ; 
lo 5 E: 
(yz|x2)q9 T V3 | (ddo), — Es (dda), — = qun, | віп Ё cos Č X 
B 1 D 1 
x SUE Б l cos Ss H 
AG 
(yz|x? — у®)» x V2 віп С ТЕ (ddo), — 5. (dam) se E un, | cos 6 sin e п — 


(ddo), | 3 (4ал), 3 (ddó), 5 3 ] 
+ — i sin ( — ( 0), — — ( ), + — ( Ji cos E cos —n + 


dd S айл 5 add а 
ck ( UTE lE ), Ses al 


mia 1 1 
(у2/32 — г?) -- 2 V5 B (ddo), — 3 (дал), — r3 un, sin сх 


5; Es 1 Rm 2 +i E 1 2 
cos cS --- --- --- -- = 
an 3 7] sin 3 H L cos ¢ cos 3 7 cos 3 7 


(х2|х2)ц d, + [(ddz), + 3 (ddó),] соз £ cos n + 2 (dd л), cos 2 € 
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1 
(xz|xz),s cos 1 [2 (ddo), + (дал), + (ddé),] cos £ cos E 9 + 
2 2 
+ = [2 (ddz), + (ddó),] cos же | + 
А 1 
+ z cos 1 [2 (ddo), + (аал), + (446),| cos E sin E 7- 
2 2 
ез [2 (Фал), + (dd9),] sin TA 
Sek cam 1 1 
(х2[322 — г?) —2V2 El (dde), — Ei (ddr); — e (ddó), | sin £ sin Ẹ x 
1 Ol 
x SE ый 
3 1 
= yt xt Юн dot: E (ddo), + 3 (ddr), + zr (dd9), | cos 6 cos 7] + 
1 
4- ex [3 (ddo), + (dd9),] cos 2 E 
(x? 2142 3) d cos Ё Е (428), + H (ddr), + C (ddà), | cos £c : + 
— — — — — --- (0) — 
M У Is 3 Т 1 3 1-49 1 55 COS 3 1] 
1 2 25 2 
+2 ФТ (ddo), + сү (Фал), + 15 (449), Сор 4- 
мыс б 1 (445) Dot E ru 
в T EET л) т ре п! 
1 2 25 — 
-2|- (ddo) + < (dda), + = (ddd), | sin =n 


(х2 — y?/3z2 — т?) B V3 ](ddo), — (d49),] (соз ё соз n — cos 2 3) 


N 


(a? — у2]322 — r?) 


idu 2 5 
V3 B (ddo), — ii (ddr) s "8 un, | cos 6 X 


1 2 ch y^ ғ 2 
x pasce ЕР. + i | сов We эсу 


(322 — г°]34® — г®) d, + * [(ddo), + 3 (44),| (2 сов $ cos n + cos 2 5) 


'8; 


1 
(322 — 72/32? — т), E (ddo), 4- = (ddr), + = (dd d cos Č X 


2 
1 2 ЖБК i 169 
— — == — sin — 
х 2 соза орал} соз d +i | 2 cos praesta 8 CH 


E 
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It should hardly be emphasized that all the above tables can be also used for the 
reduction of general overlap integrals of Bloch sums. 

The author wish to thank Dr. M. Suffezyriski for suggesting the theme of the 
present work and for valuable discussions. 


КРАТКОЕ СОДЕРЖАНИЕ 


М. Мьонсек, Применение метода тесной связи для исследования своиств симметрии 
энергетических полос в компактной гексагональной структуре. 


В работе высчитано элементы матриц энергии сперва при применении 
интегралов трицентровых а после при приближении двуцентровом. В расчетах 
было принято во внимание только первых соседей сетки. 


J REFERENCES 


Bloch F., Z. Phys., 52, 555 (1928). 

Slater J. C. and Koster G. F., Phys. Rev, 94 1498 (1954). 

Herring C., Hill A. G., Phys. Rev., 58, 132 (1940). SC 3 
Schiff B., Proc. Phys. Soc. [London], A68, 686 (1955) | 
Fletcher G. C. and Wohlfarth E. P., Phil. Mag., 42, 106 (1951). 
Fletcher G. C., Proc. Phys. Soc. [London], A65, 192 (1952). 
Suffezynski M., Acta phys. Polon., 15, 111 (1956). 


ый У чай 55. еса eh 


Vol. XVI (1957) ACTA PHYSICA POLONICA Fan 


PI-SIGMA INTERACTION IN THE FREE-ELECTRON MODEL. I. 
THEORY. 


By S. OLSZEWSKI 
Electrode Processes Laboratory, Institute of Physical Chemistry, Polish Academy of Sciences, Warsaw 


( Received April 23, 1957) 


The method of antisymmetrized molecular orbitals for the three-dimensional free- 
electron model was developed here for the case of linear conjugated molecules. In this 
model the л and с electrons are assumed to form а gas in one cylindrical box extending 
along the molecular core, and both groups of electrons are characterized according to the 
classification given originally by Nikitine and Komoss. Then the whole z — c interaction 
can be treated directly, contrary to the treatment of this problem in the Goeppert-Mayer- 
-Sklar scheme. 

The values of the electronic repulsion integrals presented by the author in an 
earlier work for the case of the so called “one-dimensional” free-electron model, were 
verified as corresponding to the integrals of the present scheme when the radius of the 
box tends to zero. 


1. Introduction 


As a further step in the development of the theory of the free-electron model, 
we may consider two problems which are quite important, if we want to undertake 
more advanced investigations of absorption spectra by means of this model. The 
problems are: 

a) what is the influence of the interaction between the more stable (steady) 
c-electrons and the z-electrons on the spectrum of the molecule; 

b) how do the transversal box dimensions affect the resulting location of energy 
states. ! e 

Before giving the mathematical scheme of our treatment of the above problems, 
let us recollect briefly the basic ideas of the л — 6 interaction calculations in the 
.LCAO MO theory. These ideas were recently discussed very keenly by Stewart (1956) 
and the present report on them is based mainly on the results of his work. 

The characteristic feature of the AS САО MO method is, that although the л 
electron states are described in this method exclusively by means of the molecular 
orbitals constructed from 2p atomic orbitals, the other wave functions of the L-shell 
are also taken into account in computing the energies corresponding to those molecular , 
orbitals. 

кт _ (369) 
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In view of this fact it is not permissible to consider the results of the LCAO 
method as related to л electrons only, for we have here ап irresolvable "intermixing" 
of z and о electron energies caused by the existance of interaction integrals between 
the л electrons and the core. Therefore in the Goeppert-Mayer-Sklar theory we cannot 
speak about the energy of a л electron system in an exact sense. It is also not possible 
in this scheme to obtain correct values for the differences of the total electronic energy 
of the molecule, owing to the inadequate determination of the core potential. 

Thus, in order to obtain a correct description of molecular electronie states the 
entire л and c electron system should be examined. The point nuclear charges with 
ls electrons might then be taken as the core. The inclusion of o orbitals into the AS 
LCAO MO scheme in order to perform л — 6 interaction calculations is of no avail, 
in such a meaning that the c electrons just play an important role in determining the 
electronic energies and eigenfunctions (through the choice of Z) of the л electrons 
themselves. 

Hence, the conclusion of Ross (1952) and Moser (1953) that the zz — о interaction, 
іп the meaning discussed there, is small enough to be neglected is not sufficiently 
established. We know only that the л — о interaction energy is to a great extent 
automatically included in the standard ASMO calculations of the electron states. 
In this scheme, we know, for instance (Stewart), that the share of the 7 — с Coulomb 
energy in the whole interaction energy in ethylene states amounts to about ten times 
that of the zz — л Coulomb and exchange interaction energies. 

So, we see that the Goeppert-Mayer-Sklar scheme — as developed by Ross and 
Moser — cannot be used at all to investigate the л — С interaction effects. Thus, 
an attempt should be made to build up a theory, in which the interaction between 
с and л electrons can be directly and fully (i. e. all the 7 — z,z — o and с — в in- 
teractions) taken into account. This corresponds to a computational scheme in which 
л and c eigenenergies and eigenfunctions are separated and completely independent 
from each other, in the sense that in the л eingenfunctions and energies there is no 
reference to c, and vice versa. 

It may be noted that each of the three types of interactions mentioned above can 
then be treated separately. 


2. Treatment of the Problem 


Nikitine and Komoss have given, in the author's opinion, a reasonable method 
of describing z and c electrons in molecules with a conjugated double-bond system 
based on the theory of the free-electron model. They have assumed, similarly to treat- 
ment of л electrons thus far in the free-electron model scheme, that the o electrons 
belonging to the single and double C—C bonds of the chain form a gas. Nevertheless, 
the both gases are characterized differently in view of their geometrical distribution 
and energetical relations. 


Let us consider, following Nikitine and dme. the solution of the Schroedinger 
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equation for free electrons іп a cylindrical box of the length L with a straightened 
chain of carbon atoms аз the axis; we assume, in addition, that the potential inside 
the box is equal to zero. The solution has then the general form: (for a non-branched 
box; the conditions that the wave functions vanish at the ends of the free-electron 
path being satisfied): 


> 9) d 
Фыл (Г) = Vi sin = z № Л (ит) cos 19 (1) 


The energy eigenvalues are: 


hè [ и? п? 
Ey, = 8m E S (2) 


(№, — normalization factor in r, д variables). 


Next it was assumed that solutions of the Л-йуре correspond to ø electrons and 
solutions of the J,-type to л electrons. Such a treatment is not only in agreement 
with the energetical classification of both electron groups, i. e. for a given n, q state 
a higher energy corresponds to the л than to the с electrons but, И has also two additio- 
nal features: 

1) — for both electron groups it is not necessary to assume a priori any spatial 
restrictions in the r direction, since these restrictions are already contained in solutions 
(1): the constants uy, and ш, can be chosen in such a way that Jọ and Л vanish 
simultaneously on the cylinder's edge, and the first zero points were taken as those 
beyond which both functions vanish. (Henceforth we shall omit the index q in the 
wave functions as equal everywhere to 1). 

2) — when u, and u, are chosen in the above manner, the eigenfunction of 
c electrons, at the point where the radial function of the z electrons (Л) attains 
a maximum, assumes approximately the same value. The radial function of o electrons 
has its maximum on the cylinder axis. Both functions, as already mentioned, vanish 
simultaneously for r — R. (R denotes the radius of the tube). 

Thus, the above solution to our problem has the following simple physical 
interpretation: the л electrons are “‘solved“ in the с electron **sea*, and the o cloud 
is most dense on the axis of C—C bonds; the л electron cloud has its density maximum 
approximately at half the distance between the axis and the side-wall of the box. 
At the place where the maximum of the л electron density occurs, the с density is, 
although only roughly, as large as that of л electrons. 

The picture described above is qualitatively in very good agreement with the 
analysis of the mutual relation between o and л densities made by Coulson et al. 
(1951). In the essential points both descriptions are, in fact, identical; hence the 
application of the Nikitine-Komoss model to л — с investigations seemed to be 
especially convenient. | | 

Most important is the fact that now we may perform the calculations for the 


many-electron problem where the л and o states are characterized independently 
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and the total mutual interaction between the electrons is directly treated. In this 
way we avoid the internal contradictions discussed in the introduction which arise 
in the calculation of z and o states when we start from assumptions involved in the 
methods based on computations with LCAO orbitals. 
We have the Hamiltonian in the form: 
Ни H 
и ues W 

where should be taken in cylindrical coordinates. The eigenfunctions in the zero- 
-order approximation foro and л electrons are those of formula (1), with indices n, 0 and 
n, l, respectively. We construct the antisymmetrized wave function for the whole 
system of both electron groups by assuming that the л and с pairs occupy, successively, 
the energy states beginning with the lowest. In the case of e. g. polyenes with М car- 
bon atoms in the chain, we have 2(N— 1) electrons in the o-group (on N—1 levels) 
and N л electrons (on N/2 levels). The electrons from the С--Н bonds are strongly 
localized and are mostly outside the л electron "sea" — hence we do not include 
them into the л — с gas of our problem. We also assume, similarly to Ross and 
Moser, that the interaction of these electrons with the electron pairs of the C—C 
bonds is not essentially significant. Thus, the wave function of our s-electron 
problem (s = 3N—2 for polyenes) will be composed of antisymmetrized products 
of the type (о, В — spin eigenfunctions): 


(91, a), (91,0 В)» (92,0 о); 4... (91.1 a). (911 p (92.1 a); es 

(Pio а)? Lë o В)? (Фә, а)? WH (914 a)? (913 В)? (P21 а)? es 

ire ект лана сан ВЕ LE 

б h ZUM [TP РЗ Sus Жез A E le Lee NUM EE A 
E pie ree Se A a Pere a О Ren oe tes (4) 

(Pio a); (P1,0 PX (P2,0 a); e (Pin el (91, В); (92,1 а) e 

(91,0 a); (91,0 В)? (фо о a); eee (913 a) (Pia В)2 (P21 а)2 Zr 


efe емо 9 ъа 6.64 a s» e À ое e 4 e e ғ 974 6 ола 2119) 929% 9.9 6 s V € 
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The energies of the states and the excitation energies can now be computed 
according to the generally known scheme; in addition, let us remember, that now 
only the л — л and с — о transitions are permissible. л — o transitions, as Nikitine 
and Komoss have shown, are forbidden. 

Now in the scheme developed above it remains only to give the formulas for the 
electronic interaction integrals. We shall have three types of езе integrals. In the 
Coulomb integrals case; for instance, si are: 


Чит = - | |; nl 6 ) Pry (г Ши 2 Pat (ғ ) Ж ш; (г ) dt dei (5) 
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and correspond for | = Г = 0 to с — о interaction, 
for = 0, Г = 1 or vice versa — to л — с interaction, 


and for l= [ = 1 — ю л— л interaction. 


We obtain three types of exchange integrals in a similar way. 

It is also possible, in our treatment, to take the configurational interaction into 
account. Then the index system of o in the electronic repulsion integrals becomes, 
in general, completely irregular. This, however, as we shall see, does not lead to any 
additional complications in the calculations. 


3. Calculation of Interaction Integrals 


Each of the Coulomb or exchange integrals may be very easily reduced — as it 
was been done previously (Olszewski 1955 — henceforth referred to as I) — to the 
calculation of integrals of the type 


JB "ій 


E Jr (шт) Л (шт) 
С == T № yi cos a z cos c z' dz dz Ji (ur) ог x 
Л (шт) Jr (ur 7) 


cos l 9” cos 1. 9 ( ed 
X Jr (up r) гг dr dr’ cos | or cos 9749492 (6) 


cos 19” cos Г 9 712 


where гь = |r — r'|, and /, / are independently 0 or 1. 
Now let us make use of the fundamental relation (Ivanenko and Sokolov 1951): 


C= [о(@Ф(®)4с= f f o(.^26(. г) o (r) ат ат - (7) 


where Ф is the electric potential of 9’ (г), and С is the Green function of а three- 
dimensional Laplace-Poisson equation, i. e. satisfying the equation 


pn le н) | (8) 


In eq. (7) о and о” are the products of the functions appearing in the C integral, in 
(6) and are dependent on г, 9, z and г’, 9°, 2” respectively. 

Next, we rewrite д (r — ғ) in cylindrical coordinates 
о : (9) 


and, further, we represent ô(z — z') and ó(r — ") in the form of the expansións: 


6 (r — r') = ô (z — z’) 
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+оо 


1 ) 
SS i x |z—2'| 10 
ô (z— 2’) әл | е d x (10) 
= d i Е. Jy (Er) Jo (kr) dk (11) 
б 
(Let us remember, that the relation ô (г — г) = Г 6(r—r’) holds). 
r 


Dividing (8) by the р? operator, we obtain finally +? 


+оо оо 

З 1 ei х |2—2'| > P 

ced ] mas | enn ana.se—n (12) 
—со 0 


Now we may introduce the expression for С into (6) and perform the integration. 
Unfortunately, it is not possible to complete the integration of C. Nevertheless we 
obtain a result which can be relatively easily treated numerically. 


We have? (for c= а): 


йай ERER 
C » cos? DEE 
a= №№ DS а wei de (OG о e 
0 о 0 


cos? (c — a) — 


— 4, (&a)) + Z -— {h (дд) — by (c) — fy (да) + L (а)} + 


1 The absolute value |= — z'| in the Fourier expansion of ó(z — 2’) ensures that the sign of the 
exponential coefficient does not change and thus makes possible the integration. 


2 [t is a fact of essential importance that when the full ô (r—r’) expansion is used we may interchange 
the order of integration over z, z' with the integration over k (and over x) as a result of the presence 
of the product Jo (kr) Jo (kr"); then the double integral over х, kin С is uniformely convergent in the 
whole z, z' region. However in the AS FEMO method, previously developed in I, we have also made calcu- 
lations with the three-dimensional wave functions of the FE model, and of course, with the corresponding 
Green function for the interaction potential, but it was assumed next, that the potential acts in one 
direction, i. e. along the free-electron path. Then the Green function takes the form: 


Сб sd Seog ast [oe kak |49 | 
o 4л|@—4] 8m J x? + k? = (13) 


Here, however, С”, is divergent for z = 27 and the interchange of the order of integration over k, x 
. with that over 2, 2' is not permissible. However, see further the next section. 


3 The sign of C corresponds to the absolute values taken for the interaction integrals. 
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c? 
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where we have introduced the following notation*?: 


© = (r? 4-r'2—2rr' cos y)” 
l v4-2m 
— @а 


=}, элк теа 


— Bessel function of the second kind 
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Y »+2т +1 
— фа 


Е, (ба) = м Теа Т — Struve function of the second kind 
1 2 1 
М (ба) - m S m x n(—1)" E Qa 
H m 


(tor large a, M (a) tends rapidly to — а. Lai (ба) — 1] үн 


For the Bessel function product in { e (and for the cosine productin[ ] in Eq. 
(6)), we take upper or lower row for Coulomb or exchange integrals respectively. 
D is a constant equal to 


2x when [= l'=0 (o — o interaction) (18a) 
abun ue os егіріреісіелі (18b) 
Әл when 1 = Г = 1 _(л — л interaction) (18с) 


In the computation of the above integrals the ойыр of the wave functions at the 
ends of the free-electron path was taken into account. In addition, we note that when 
the configurational interaction is included in our scheme, there may appear in the 
matrix elements integrals of a type different from those calculated above. They have, 
under the integral sign, the form of the following expression corresponding to the 
product of "transversal" factors of the wave functions of our problem: 


Jing) Dos (nr) Je (nr!) Je (ptt) nos 10 cos PO! cng FO Ros ee 


where J, I’ can be 0 or 1 independently one of another. Such integrals, as may readily | 
be seen, vanish for lÆ Г as а result of the integration over 9, 9’. - Thus only ішін | 
Ste E actic 
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the molecule and then the experimental values on the C—C distances are taken into 
account, the second one can be estimated either on the basis of the considerations of 
Nikitine and Komoss, or from formulae given by Kuhn and Huber (1953). In the latter 
case К is connected with the effective nuclear charge of the carbon atoms, being (R) 
inversely proportional to Z. However, for larger molecules, the criteria of Nikitine 
and Komoss lead to contradictions. Besides, in their method the values confining R 
change considerably with a change in the chain length; this is especially so in the case 
of polyenes. 

In our treatment of the problem it seems most reasonable to take Z equal to 
four. For such value of Z the radius of the box would be confined, according to Kuhn 
and Huber, to lie between 1.15À and 1.35À, assuming that the side-wall is placed on 
the boundary of z electron cloud. This value of R is in very good agreement with 
that obtained by Nikitine and Komoss for the first homologues of symmetrical cyanines. 


4. The case of R tending to zero 


The result which we obtain for the interaction integrals when the radius (or the 
ratio of the radius to the box length) tends to zero is significant, since in this case 
C attains the values obtained in the "one-dimensional" AS FEMO scheme. 

Previously (see I) the components of the interaction integrals were given in the form 

Sab 


2 
E 4; А) ШЕ ар Z COS а, 2' * SEN dz dz (20) 


оо 


which, however, was erroneous (see footnote 2), since the integrals of the (20)-type 
are divergent. Thus we have a case in which the un: of the limit, i. e. 


ШІ e (г) о’ (Г): —— Tux EE (21) 


and hence, in view of normalization of g іп г, 9 я the integral :(20)®7, is 
not equal to- the limit of the integral ч 


im (ament. = incu Ms i (22) 


In the latter case we obtain the results given in Н (See Appendix for the proof)®. 


6 We have assumed here, that the parts of o depending on z have the more general form of the type 


1 1 
A, sin а, 2 (A, is the normalization factor in 2). Then in C (see (6)) we have instead of та the factor —Ap A 


7 Таш much indebted to Dr. Н. Labhart for drawing my attention to the divergence of the inte- 
gral (20). 
8 We obtain this result when for D we х the value 27 in the C calculations of the present work. 
This corresponds to an equal distribution of the electron density about the box axis, hence to the as- 
sumptions of simple FE model. 
- # 
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'The conclusions presented recently have interesting consequences from the point 
of view of classical electrostatics. This concerns to the problem of the interaction 
energy of charges in a very thin tube. 

If we assume a priori, that the interaction potential in such a filament is of the 
type 

1 


8—2 


(23) 


then after integration we obtain for the interaction energy a quantity tending 
logarithmically to infinity. On the other hand, when the problem is treated from 
the beginning as a three-dimensional one, and when with the result for the interaction 
energy we next pass to the limit for the cross-section tending to zero, we then obtain 
a finite value. For instance, in the case of two electrons moving freely in a filament 
e? 

of length L, the interaction energy attains 1:66 the value of integral given in I in 
the case of c — а = 0. 

This second interpretation of the calculation of the interaction energy seems to 


have stronger physical grounds. 


Appendix 


It is essential to give the proof that as К — 0 the expressions obtained for the 
interaction integrals in the scheme developed above, go over into expressions ob- 
tained in the AS FEMO treatment (in I). Now let us perform those calculations for 
the simplest example: C in the case c — a — 0. The proof for other cases in quite 
similar. 

We majorize by unity the Bessel functions appearing in the components of the 
C integral with minus signs (see Eq. (16) of this work). We shall show further that 
these integrals tend to zero for R — 0. ` 

In the first negative term we replace © by the series: 


oo 


D Ж Са” (cos x) (5) ‚(5 < n) (24) 


where C; аге Gegenbauer’ s coefficients of the expansion of ©. After integration 
over г, г! we obtain for the first negative term the expression: 


Т л 
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0 
Integration over у gives quantities not larger than 2. The normalization product is ` 
proportional to R74, hence the whole expression is proportional to R. Thus it vanishes | 
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when К tends to zero. Similarly the second negative term (with J equal unity) gives 
after integration (here, for a = 0, іп M (wa) expansion we have m = n):? 


S » Dis —1])/ xo (mV 1 2m 
9 (са Gah (t) То ** (26) 


(А, В denote the constant factors in front of the integrals). Again taking the norma- 
lization factors into account, we find that the last expression is at least of the second 
order іп К, and therefore also vanishes when R -> 0. 

Thus there remains the first positive term in which the integration over г, r’ with 
respect to the corresponding normalization coefficients of J functions gives simply 


2 
о $ р 5 E : д € Rots uc о 
unity. Taking all the factors into consideration we obtain C equal to т. This is identical 


with the value of C obtained for the same case in I (after dividing the latter by L?, 
since in the previous definition of C the normalization factors of the functions de- 
pending on z were neglected). 


КРАТКОЕ СОДЕРЖАНИЕ 


C. Ольшевски, Воздействие m—6 в модели свободного электрона. 


Развито метод антисиметризованных молекулярных орбит в тримерной мо- 
дели свободного электрона для линейных сопряженных молекул. Электроны ли б 
создают газ в одном цилиндре растянутым вдоль цепи молекулы, а обе электрон- 
ные группы схарактеризированы согласно с классификацией поданной Ники- 
тиным и Комоссом. Тогда полное воздействие л и с может рассматриваться 
прямо, противоположно к подходу к этой проблеме в схеме Гепперт--Маер- 
Склар. Результаты на интегралы электронного воздействия в случае т. наз. 
‚ одномерной модели свободного электрона, поданные автором в более ранной 
работе, были верификованы ‘как отвечающие интегралом полученным в сейчас 
рассматриваемой модели для случая когда радиус цилиндра приближается 
к нулю. | 
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9 Тһе term by term integration is permissible here in view of the uniform convergence of the 
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Ап investigation was made of nuclear magnetic resonance in a stream of flowing 
liquid. The dependence of the amplitude of the absorption lines on the rate of flow was 
investigated. It was established that the flowing liquid method allows a simultaneous 
determination of the saturation factor and the longitudiual relaxation time; in this con- 
nexion, the relaxation time measurement may be made for any (not too small) degree of 
saturation. It was also shown that the limit of applicability of the method can, in comparison 
with the limit of applicability given by Suryan, be shifted considerably in the direction 
of long relaxation times. 


Introduction 


Suryan (1951) made measurements of the amplitude of absorption lines of nuclear 
magnetic resonance in flowing liquid as a function of the rate of flow. If the following 
conditions are fulfilled: 

a) the liquid, before entering the resonance coil, is in the magnetic field long enough 
for stationary magnetization M, = xy Ну to be attained, and b) the saturation of the 
resonance line is sufficiently strongt, it should be expected that there will be a 
linear increase of the line amplitude with an increase in the rate of flow 


D 

where A is the observed amplitude of the resonance line for а rate of flow v, A, is the 
amplitude of the saturated line in the liquid at rest, and / is the length of the т. f. 
coil. Suryan obtained experimentally for a solution of N/1000 FeCl, in water a linear 
dependence of А on v, and from the slope of the straight line he determined Ti 
— 0.06 sec. í 
——— 

1 Suryan gives the: condition 1 y? H? T, T,* >> 1. Since the amplitude of modulation of the 
magnetic field H,, in his measurements is greater than the line width, this condition should have the 
form у H?T,/H,, >> 1. 


: (381) 
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The flowing liquid method allows a considerable improvement of the signal-to- 
-noise ratio, which is of great importance in the case of small signals. For this reason, 
Denis et al. (1952), Bloom and Shoolery (1952), and Mitchell and Phillips (1956) 
used this method. 

Gaussen (1955) used another method for measuring the time Ту by means of 
flowing liquid. The liquid, without previous magnetization, flowed through a resonance 
coil. After a sudden stopping of the flow it was observed that resonance lines arise 
as a result of the magnetization of the liquid which stopped in the coil. The recording 
of this process on film permitted the determination of T}. 

The flowing liquid method was used by Sherman (1954) to investigate nuclear 
magnetic resonance for separate regions of excitation and detection. 

In this work, a more detailed investigation was made of the dependence of the 
amplitude of the resonance line on the rate of flow of the liquid in the stationary 
state. It was established that the flowing liquid method allows simultaneous determi- 
nation of the saturation factor, and, in this connexion, the measurement of the rela- 
xation time can be made for any value of this factor. It was also shown that the limit 
of applicability of this method can, in comparison with the limit given by Suryan, 
be shifted considerably in the direction of long relaxation times. 


Simplified Theory 


The process of the saturation of samples by a high frequency field is discernable 
in the decrease of amplitude А of the absorption line and is described by the following 
equation : 


t 


A zd: = (24222224) е PU , ` (2) 


where A, is the initial amplitude of the line (at the moment the r. f. is applied), A, is 
the amplitude of the line after saturation, and z = А,/ А, is the saturation factor 
dependent on the amplitude of the г. f. field H}, on the relaxation time Т; and on the 
amplitude of modulation of the magnetic field. 

We denote 


Е UTD 9) 
The liquid flowing with a speed v flows through a coil of length [їп a time t = /,. 
We assume that before entering the resonance coil the liquid is in the magnetic field 
long enough for stationary magnetization M, = x, Mto take place. The mean value - 
of S for the region of the resonance coil is 


I ЕТЕ oT (125 iE. 
5 = ТЭ fe »T,z dy — e, "n ens ). (4) 
0 Ф 
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There should be taken into account the velocity distribution of the liquid as a function 
of the distance r from the tube axis 


ЕА o 


where R is the tube radius and v is the mean velocity of the liquid. The averaging 
over the cross section of the tube leads to an integral which cannot be evaluated 
analytically. It is sufficient, however, to restrict oneself to the examination of the 
limiting cases: 


l 


v Tiz 


> 1 (6) 


This condition corresponds to а slow flow of the liquid or to strong saturation. 


The neglecting in Eq. (4) of the term exp (—//о Т, z) leads to the formula 


ET (7) 


l 
EC (8) 
This condition corresponds to a rapid flow or to weak saturation. Expanding 
exp (—1/v T, z) into a series and neglecting terms of the order higher than опе, 
we obtain 


Ж ee 


(9) 


5 
2 


Тһе fulfilling of condition (8) may be іп doubt, since the averaging over the 
tube cross-section gives a divergent integral. However, it is readily shown that this 
does not create any basic difficulty. For, integrating over the tube radius in the limits 
of from 0 to 0.99R, and therefore, neglecting barely 0.2% of the volume of the flowing 
liquid, we obtain = ^ = = z , which causes very little difficulty in 

No А ФТ! = 
satisfying condition (8). 

The solution given leads to the conclusion that for a low rate of flow we have 
a linear dependence of s on v while for sufficiently large rates of flow s ееазез to depend 
on v and attains a constant value s,. The greater the coefficient z, and hence the weaker 
the saturation of the liquid in the stationary state the smaller the value v for which 
the value s, is attained. The line amplitude corresponding to the value of 5,, as may 
be seen from Eq. (9), is equal to Ay. Therefore, from the determined values of Ag 


d 
7 
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and А, we can calculate the saturation factor 
== 5 (10) 


and then, from the slope of the curve of А as а function of v we can calculate the 


relaxation time T, 
A — A; 1 


Ийт DEET (11) 


А certain systematic error may be introduced by the fringing effect in the coil, 
owing to which the length is not precisely defined. This error, however, does not 
exceed a few percent, and in making the relative measurements, it was completely 
eliminated. 


Equipment and Experimental Procedure 


The measurements were made for liquid containing protons by means of a bridge 
arrangement at a frequency of 28 Mc/sec in a magnetic field of 6,580 Oe. The electro- 
magnet used had a pole face diameter of 140 mm and a gap width of 30 mm. An 
energizing current of 5.4 А was supplied by an accumulator battery. By means of 
three differential screws which allowed the pole faces to be adjusted precisely parallel 
to one another a change in the magnetic field of AH — 0.02 Oe was obtained in a region 
of the sample 15 mm in length and 7 mm in diameter. The block diagram of the 
equipment is shown in Fig. 1. An г. f. bridge of the double — T type given by Grivet 


| Ho * Hin COS W mt 


rf 


receiver 


22615 


photocamera 


Fig. 1. Block diagram of the equipment 


et al. (1951) was used. The liquid flowed through a glass tube, the inner diameter 
of which was accurately measured. On the tube was wound the r. f. coil. Before enter- 
ing the coil, the liquid flowed through a helical tube placed between the pole pieces 
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of the electromagnet. This allowed the stationary value of magnetization to be attained 
(Fig. 2). The length of the path of the liquid in the magnetic field was 360 cm. or, 
with a flow of 40 cm/sec, the mean time spent by the liquid in the magnetic field 
was 9 sec. This time, which in the case of water is 3.5 times longer than the relaxation 
time, was sufficient, in practice, for the full value of M, (М = 0.97 My) to be attained. 
In the case of a liquid with a shorter relaxation time, the magnetization was still closer 


| intet 


torf bridge N 


outret 


Fig. 2. Path of liquid between pole faces of the electromagnet 


to the value Mọ. The liquid flowed under a pressure of about 0.25 atm from a reser- 
voir placed under the ceiling of the laboratory. In order to facilitate the change of the 
rate of flow a set of 20 capillary tubes with suitably chosen diameters was used. These 
were placed at the outlet of the liquid stream. The capillary tubes were calibrated 
by accurately measuring the rate of flow for each. Such set of capillary tubes consider- 
ably facilitated the experimental procedure, since the measurements could be repeated 
several times in short intervals of time without the necessity of making each time 


Fig. 3. Photograph of resonance line of distilled water for different rates of flow a) v= 0.5 cm/sec 


f b) v = 40 cm/sec. 
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a time-consuming measurement of the rate of flow. For various rates of flow the 
absorption lines were photographed from the oscillograph screen. In Fig. 3 are shown 
sample photographs of absorption lines in destilled water for two different rates of 
flow. In order to avoid the influence of any eventual detuning of the bridge, high- 
frequency generator, or receiver, every second measurement was made for a standard 
capillary tube and the line heights as measured from the photographs were referred 
to the line heights for this capillary. 


Results of the Measurements 


Several series of measurements were made for ordinary distilled water and for 
ап n/1000 solution of MnSO, in water. For both liquids measurements were made 
for two different values of the high-frequency field, and hence, for different saturation 
factors. . 


Figs. 4 and 5 show curves of the line amplitudes as a function of the rate of flow 
of the liquids. The shape of the experimental curves fully confirm the theoretical 


Distilled water 


Resonance line amplitude 


2 (0 20 30 40 7 
: Flow velocity (cm/sec) 


Fig. 4. Plot of amplitude of resonance line vs rate of flow for distilled water 


predictions. From the curves drawn on the basis of Eqs. (10) and (11) the saturation 


factors and the longitudinal relaxation times were calculated. The following values 
were obtained: | 


Distilled water | 
H, lower = 15.8 Т, = (2.60 + 0.20) вес _ 
H, higher Z = 1/98 T, = (2.77 + 0.25) sec 
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1/1000 n MnSO, solution 


H, lower Z = 12.2 Т, = (0.246 - 0.20) sec 
Hi higher ZI 113-0 Т, = (0.257 - 0.20) sec 


The errors listed were estimated on the basis of the accuracy with which the values 
of Ay, A, and the slope of the linear portions of the curves may be read. 

The results of this work show that the flowing liquid method can be used for 
measuring the relaxation time also when the saturation is not large, since the line 


MnS04 solution 


Resonance line amplitude 


SI 
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Fig. 5. Plot of amplitude of resonance line vs rate of flow for n/1000 solution of M, SO, in water 


amplitude measurement for high rates of flow permits the direct determination of 
the saturation factor. As may be seen from the example of distilled water, the method 
described can be used for measuring long relaxation times of the order of seconds. 

We would like to thank Professor H. Niewodniczariski for making this work 
possible in his laboratory and K. Krynicki and G. Zapalski for valuable aid in making 
the measurements. 


КРАТКОЕ СОДЕРЖАНИЕ 


А. Грынкевич и T. Валюга, Ядерный магнитный резонанс в текущей жидкости 


Проведено исследование ядерного магнитного резонанса в текущей жидкости. 
Определено зависимость амплитуды линии абсорбции от скорости течения. До- 
казано что метод текущей жидкости позволяет одновременно получить коэф- 
фициент насыщения и продольное время релаксации и поэтому измерение вре- 
мени релаксации можно исполнить при любой (не слишком малой) степени насы- 
щения. Показано тоже что предел измеряемых этим методом времен релаксации 
по сравнении с пределом поданым в работе Сурияна можно значительно расши- 
рить в направлении длинных времен. 
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The authors investigated the operating conditions of a dielectric resonant amplifier 
and the dielectric properties of (Ba—Sr)TiO, ferroelectric condensers employed in the 
amplifier. Curves of the relative permittivity variation versus the ОС field at various 
temperatures were obtained. The following characteristics of the dielectric amplifier under 
investigation were measured: output voltage versus signal voltage and signal frequency, 
and voltage amplification coefficient versus temperature. The largest voltage amplification 
value was obtained within the Curie point region. The temperature being constant, this 
amplification value remained practically constant over the entire audio frequency range. 
The measurements proved that it is possible to determine the optimal operating conditions 
of the dielectric amplifier from the graphs showing the dependence of the permittivity 
of the ferroelectric material on the DC field and on the temperature. 


Introduction 


During recent years a number of papers (Vincent, 1951, Pipes, 1952, Urkowitz, 
1952, Penney et al. 1953, Silverstein, 1954, Mason and Wick, 1954) have been 
published on the use of ferroelectric materials as nonlinear elements in dielectric 
amplifiers. Since the results hitherto published are insufficient for determining the 
optimum operating conditions of an amplifier, the authors submitted the dielectric 
properties of ferroelectric materials employed in amplifiers, ‘and the operating 
characteristics of dielectric resonance amplifiers, to a detailed investigation. 


1. Preparation of ferroelectric material 


The ferroelectrics investigated were obtained from solid (Ba—Sr) TiO, solutions 
by known methods (Trzebiatowski et al. 1952). The authors prepared a number of 
series of ferroelectrics varying in Sr content, starting from chemically pure BaCO}, 
SrCO, and TiO,, which after being mixed, were sintered for З hours at 1300°С. The 
(Ba—Sr) TiO, thus obtained was ground and subsequently compressed into disks 


(389) 
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12 mm in diameter under pressures of 5 — 8 t/cm?. The disks were then sintered 
once more for 2 hours at 1200?C and 3 hours at 1350?C consecutively, yielding material 
of low porosity. After the disks had been polished down to a thickness of 0.15— 
0.20 mm, the surfaces were covered with silver electrodes by evaporating in vacuo. 


2. Measurements of field and temperature dependence of the dielectric permittivity 
in (Ba—Sr) TiO, 


The dieletric permittivity in polycrystalline (Ba—Sr)TiOg falls with growing field 
bias (Piekara and Pajak, 1953). Using the circuit shown in Fig. 1 it could be proved 
that the conditions prevailing during operation of a dielectric amplifier are the 
following: as the signal voltage increases, the dielectric permittivity of (Ba—Sr)TiOg 
decreases, and conversely, during the period of decreasing of the signal voltage the 
permittivity increases. Low frequency signals from an audio generator were supplied 


VALVE 
SWITCH 


HUE 
GENERATOR 


AUDIO FREQ. 
GENERATOR 
Fig. 1 Dielectric resonance amplifier circuit. 


to the dielectric amplifier, changing the capacitances of polarized ferroelectric 
condensers in the resonant amplifier circuit and thus causing the latter to detune. 
The amplifier signals and the low-frequency generator signals were supplied to the 
oscillograph over a valve switch. Thus, the output and input signals were observed 
simultaneously. In order to ascertain whether the changes in the permittivity of 
(Ba—Sr)TiO, caused by momentary changes in the signal voltage supplied to the 
amplifier are the same as those caused by a DC voltage, the permittivity was measured 
as a function of DC voltages at various temperatures. The capacitance of the ferro- 
electric elements of the amplifier was measured by the bridge method ata теч 
of 2 Кс/з and a measuring field of 10 V/cm. 
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The results are plotted in Figs. 2 and 3. The curves show the effect of the applied 

әкетеді Sie 

field on the relative permittivity ZË at various fixed temperatures, and the permitti- 
€ 


уйу — temperature dependence at various values of the field intensity. It will be 
noticed that the field causes the greatest changes in the vicinity of the Curie point. 
It will be shown, that, indeed, in order to obtain maximum amplification, the Curie 
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Fig. 2 Relative dielectric permittivity of (Ba—Sr)TiO; versus biasing field at various const ant temperatures 
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Fig.3 Dielectric permittivity of (Ba—Sr)TiO, as function of the temperature at various biesing fields. 
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temperature of the ferroelectric condensers employed should not differ considerably 
from the operating temperature. 

The curves plotted in Fig. 2 show identical slope in the vicinity of the Curie 
point (30°С). Within this range of temperature, amplification is independent of 
temperature changes in the operating conditions of the amplifier. The maximum steep- 
ness near the Curie point (25°C — 36°C) is about 0.07 per 1 kV/cm, and falls to 0.06 
per 1 kV/cm between 20°С — 25?C and 36°С — 45°С. The ratio of these two values 
is 0.85, corresponding to 1.42 dB relative amplification units, which practically lies 
below the threshold of perceptibility in an audio amplifier. Operating conditions 
of a dielectric amplifier are considerably worse at temperatures more distant from 
the Curie point. For example, a change in temperature from 15°С to 20°С corresponds 
to a loss of the relative amplification of 4.1 dB, which is sufficient to be perceptible. 


3. Operating characteristics of the (Ва—5т) TiO, dielectric amplifier 


Тһе resonance amplifier was constructed according to Fig. 1. The ferroelectric 

condensers C, in the resonant circuit were polarized over the resistance R by employing 
d SÉ Ge 1 

a DC field of 3 to 5 kV/cm; within this range и = — Be is constant. The frequency 
€ 


of 580 kc/s which was applied to the amplifier was chosen so as to provide for operation 


Fig. 4 Output voltage U, as function of the signal voltage U . 


on the resonance curve slope of the LC circuit. Audio signals of a frequency of 0 — 20 
kc/s were fed to the coupling condenser C, and the HF filter Г,. The LF and HF vol- 
tages were measured with a valve voltmeter and a scaled oscillograph, respectively. 
The following characteristics of the dielectric amplifier under investigation were 
investigated : 


Investigation of Operating Conditions of a Dielectric Amplifier 59 


Co 


1. Output voltage versus signal voltage and signal frequency, 

2. Voltage amplification coefficient versus temperature. 

Measurements of output voltage versus signal voltage (Fig. 4) with the resonant 
amplifier operating in the vicinity of the Curie point (30°С) yielded the maximum value 


of the voltage amplification coefficient Ё = Uo = 1.3. 


0, 
It was found that the output voltage is independent of the signal frequency within 
the audio frequency range investigated (continuous curve in Fig. 5). From measure- 
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Fig. 5 Output voltage U) and power gain coefficient 7 as function of the signal frequency. 


Fig. 6 Voltage amplification coefficient k of the amplifier shown in Fig. 1, asa function of the temperature. 
Dashed curve: the coefficient k calculated from eq. (3), with the value of и obtained from the slope of 
the static characteristics shown іп Fig. 2. 


ments of the signal voltage, output voltage and also input and output impedances, the 
power gain coefficient 7 was calculated. This power gain coefficient was plotted 
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versus frequency (dashes in Fig. 5). The maximum power gain within the LF range 
amounted to abount 40 dB. 

The temperature dependence of the voltage amplification coefficient k is shown 
in Fig. 6 (continuous curve). It will be noticed that maximum voltage amplification 
occurs in the vicinity of the Curie point (30°C). Variations in temperature amounting 
to 5°C — 10°C when operating near the Curie point produce changes in the amplifi- 
cation coefficient so restricted as to be practically negligible. However, a considerable 
decrease of k appears when the amplifier operates at temperatures remote from the 
Curie point. The authors explain this decrease in the amplification coefficient by 
the fact that the largest transient variations of the permittivity of the ferroelectric 
material appear within the region of the Curie point. Hence, the Curie temperature 
of the ferroelectric material employed should lie in the proximity of the operating 
temperature of the amplifier. Similarly, the largest variations of the permittivity of 
the ferroelectric material are caused by a biasing DC field near the Curie point (Fig. 2 
and 3). 

The dashed curve in Fig. 6 shows the temperature dependence of the voltage 
amplification coefficient & calculated according to the theory of Н. Urkowitz (1952), 
using и obtained from the static measurements (Fig. 2 and 3). The change of the 
HF amplitude voltage U, as function of the signal amplitude U (eq. 39 in his paper), 
using our notations, is given by: 


0, = 0.17 Об LC (“| 0, (1) 
where / is the constant amplitude of the HF current, R — the resistance of the circuit, 
О, L, C' — resonance circuit constants, и = E um s — the thickness of each 

5 


of the ferroelectric condensers; С’ denotes the capacity seen by the signal. Denoting 
the capacity of the resonance circuit seen by the HF generator by C, we obtain 
C' — 4C. The voltage amplification coefficient is now: 


D и 
k= g —04IRo*LC (4) (2) 


In this equation w? L C = 1. For optimal conditions I R = U Үз), where U is 
the amplitude of the HF voltage. Denoting = = E, we have for the voltage amplifi- 
5 


cation coefficient: 


k=0.96 E uQ (3) 
where E, denotes the HF field strength in kV/cm, Q — the figure of merit of the 
resonant circuit and и = = га in cm/kV. The values of и for various temperatures 


5 


аге calculated from the curves given in Fig. 2 with DC fields of 3 to 5 kV/cm. For the 
amplifier investigated E, amounted to 0.75 kV/cm and О = 22.8. On comparing both 
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curves it will be seen that the variations of the permittivity of the ferroelectric material 
induced by the АС field of audio frequency, with а DC field superimposed on it, are 
almost the same as the variations of the reversible permittivity caused by a DC field 
(Fig. 2). As both curves in Fig. 6 are approximatelly parallel, it is thus possible to 
determine the optimal operating conditions of the amplifier from the DC-field-permit- 
tivity and the temperature-permittivity dependence of the ferroelectric material. 
As follows from the graphs in Fig. 6 the error in determining the voltage amplification 
coefficient by this method does not exceed about 8% 


Conclusions 


The dielectric resonance amplifier yields the largest voltage amplification 
(however, onby slightly exceeding 1) when operating within the region about the 
Curie point. The temperature being constant, the amplification remains practically 
constant over the entire audio frequency range. 

The optimal operating conditions for the dielectric amplifier may be determined 
from the curves showing the effect of the DC-field and of the temperature on the 
permittivity of the ferroelectric material. 

The authors wish to express their thanks to Professor A. Piekara for directing 
the present investigation and for his valuable suggestions and advice. 


КРАТКОЕ СОДЕРЖАНИЕ 


T. Краевский, E. Петржаки A. Сувальский, Исследование по рабочим условиям 
дизлектрического резонансного усилителя. 


Авторы провели исследования по рабочим условиям диэлектрического резо- 
нансного усилителя и диэлектрическим свойствам сегнетоэлектрических конден- 
саторов содержащих (Ba-Sr)TiOs, применяющихся в усилителе. Определены кри- 
вые зависимости Деје от напряженности постоянного поля для различных темпе- 
ратур. Для исследуемого усилителя определялись: зависимость напряжения на 
входе от напряжения и частоты сигнала и температурная зависимость коэффи- 
циента усиления напряжения. Максималньое усиление напряжения появлялось 
в области точки Кюри; для данной температуры оно оказалось практически 
постоянной величиной во всём акустическом диапазоне. Результаты измерений 
производившихся авторами приводят к заключению, что оптимальные рабочие 
условия диэлектрического усилителя определимы на основании хода зависимости 
диэлектрической проницаемости сегнетоэлектрического материала от напряжен- 
ности постоянного электрического поля и от температуры. 
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